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SUMMARY 


A method for computing three-dimensional flow in supersonic inlets 
is described. An approximate set of goveiming equations is given for viscous 
flows which have a primary flow direction. The governing equations are 
written in general orthogonal coordinates. These equations are modified in 
the subsonic region of the flow to prevent the phenomenon of branching. 
Results are presented for the two sample cases, a Mach number equals 2.5 
flow in a square duct, and a Mach number equals 3.0 flow in a research jet 
engine inlet. In the latter case the computed results are compared with 
the experimental data. A users* manual is included. 



INTRODUCTION 


The design of an inlet for a supersonic aircraft is a difficult task in 
view of the wide operating range over which good performance is desired. 

On supersonic airplanes the inlet design is compromised by the requirement 
that the plane operate at all speeds from zero to the supersonic design point, 
including the low speed, high angle of attack regime encountered during take- 
off and landing. In addition most supersonic aircraft must have a subsonic 
cruise capability where fuel economy often becomes a significant factor. As 
a result of these varied and sometimes conflicting requirements, the design 
of a supersonic inlet is a compromise. The optimum geometry at one operating 
point must be modified to be acceptable at other operating conditions. Even 
variable geometry inlets, which are required at high speeds, are limited by 
the system penalties associated with their added weight. 

There are two basic aspects of inlets which greatly influence the net 
thrust of a propulsion system: (1) pressure recovery, and (2) distortion. 

The net, or useful, thrust of an engine is the gross thrust out of the 
nozzle minus the inlet ram drag. At high Mach numbers the inlet ram drag is 
large and the net thrust thus becomes the small (relatively) difference 
between two large numbers. At a gross-to-net thrust ratio of 5, a 1% loss 
in gross thrust, therefore becomes a 5% loss in net thrust. Since inlet 
pressure recovery factors directly into the gross thrust, the pressure 
recovery is especially important. Distortion is also important, but in 
another way. When the level of distortion at the engine-inlet interface is 
smaller than the engine distortion tolerance, the effects of distortion are 
seen only in the pressure recovery. However, when the inlet distortion becomes 
higher than the engine tolerance level, the engine can no longer operate at 
that flight condition. Consequently, either the power setting must be reduced 
or the angle of attack or yaw must be reduced. Thus pressure recovery in the 
inlet is related to overall losses in the propulsion system while distortion 
is related to the necessity of operating at other than the desired flight 
condition. 

Viscous effects in the inlet are important in determining both the pres- 
sure recovery and the distortion. Three categories of viscous effects in the 
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inlet are of interest: (1) the simple boundary layer growth along walls, 

(2) bleed of the boundary layer air, and (3) shock-boundary layer interaction. 
The simple boundary layer growth along the walls of the inlet (ramp, cowl 
and sidewalls) contributes to the pressure loss in the system (as well as 
modifies the flow within the duct) by the blockage caused by the displacement 
effect of the boundary layers. Boundary layer bleed is used to control shock 
wave boundary layer interactions such that losses are kept to an acceptable 
level. Since the bleed air comes form the boundary layer, viscous effects 
are important in determining how the flow field responds to the bleed. 

The shock wave boundary layer interactions within an inlet can cause both 
pressure losses and distortion. There are typically three types of shock 
wave boundary layer interactions within rectangular inlets. The first type, 
called incident shock interactions, occurs on the cowl and ramp surfaces. 

These incident interactions are characterized by an incident shock hitting 
the wall boundary layer and generating a reflected shock. Under certain 
conditions, boundary layer spearation can occur thus giving rise to a stream- 
wise recirculation zone and large loss levels. The second type, which could be 
termed a glancing shock interaction, takes place along the sidewalls where 
the shock wave travels across the sidewall. The boundary layer in this type 
of interaction does not usually lead to a streamwise recirculation but rather 
contributes to the cross flow. The third type of interaction takes place in 
the corners. In this case, the incident and glancing interactions coalesce 
along with the merged floor and sidewall boundary layers. 

The foregoing considerations suggest that inlet design technology would 
benefit from a detailed and accurate flow field calculation procedure which 
would include the effects of the boundary layer. In addition, since the 
glancing shock wave boundary layer interactions along the sidewalls as well 
as corner effects give rise to highly three-dimensional flow fields, the 
calculation procedure should be three-dimensional. Calculation procedures 
which are either two-dimensional or do not account for boundary layer effects 
have a more limited applicability. The matching of the three-dimensional 
inviscid ’core’ flow with a three-dimensional boundary layer is difficult to 
do in a consistent manner particularly with a rotational core region. In 
addition the glancing interactions along the sidewalls and the flow in the 
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corners do not fall within the concept of conventional three-dimensional 
boundary layer theory. 

Although a full three-dimensional Navier-Stokes procedure would provide 
the necessary generality to predict the flow in an inlet, the required 
computer time and storage (Ref. 1) Indicate that such a procedure should be 
used only if no suitable alternative exists. An optimum analysis would 
possess the general three-dimensional viscous nature of the Navier-Stokes 
equations, but would take advantage of realistic physical approximations to 
limit the computer running time and storage requirements associated with the 
solution of the complete three-dimensional Navier-Stokes equations. In this 
study, such a simplifying approximation was used to reduce the complete 
three-dimensional Navier-Stokes equations to a form which could be treated 
as an initial-boundary value problem and solved by forward marching in space. 
The assumption made is that a ’primary flow’ direction exists and that diffu- 
sion arising from the rate of change of flow in this ’primary flow’ direction 
can be neglected. In this manner a set of steady state equations is produced 
for entirely supersonic flows which can be solved by an efficient spatial 
marching procedure. In any embedded subsonic regions, such as at no-slip 
walls, further approximations are required to allow solution by spatial 
marching . 

The remainder of this report will discuss the application of a three- 
dimensional forward marching procedure to calculate the flow fields in rectan- 
gular supersonic inlets, including viscous effects. The report is divided into 
five major sections. In the first section the basic set of coupled nonlinear 
partial differential equations which are derived from the Navier-Stokes 
equations will be presented. The second section will discuss the phenomenon 
of branching, which can occur in spatial marching procedures when embedded 
subsonic regions are present in the flow, and show how that problem is treated 
in the present study. In the third section the linearizations and numerical 
technique used to solve the governing equations will be discussed. This will 
include the application of boundary conditions. The fourth section will 
discuss the results of sample cases, and the fifth section will be a users’ 
guide for the computer code used to obtain the solutions. 
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LIST OF SYMBOLS 


Symbols 


a 

A 

b 

B 


c 

D 

h 

h 

o 

V 


M 

n 


m 


N 

P 

S 


u 

/-O 

u 

u 

T 

V 

X 



Parameter in Roberts’ transformation, Eq. (A-1) 

Matrix containing marching direction coefficients 
Parameter in Roberts' transformation, Eq. (A-1) 

Constant in logarithmic law of the wall, nominal value 5.0 
Parameter in Roberts’ transformation, Eq. (A-1) 

Matrix containing cross plane coefficients 

Damping factor 

Metric coefficient 

Stagnation enthalpy 

Mixing length 

Reference length 

Mach number 

Unit vector in direction of Mach line 
Niimber of grid points 

Static pressure; normal probability function 

Column vector containing source terms 

x-direction (or primary direction) velocity 

Velocity parallel to a wall 

Wall shear velocity 

y-direction velocity 

z-direction velocity 

Streamwise or primary direction in orthogonal coordinates 
Transverse direction in orthogonal coordinates 
Distance from a wall 

Nondimensional distance from a wall (Eq. 12) 

Constant in damping factor (Eq. 11) 

Spanwise direction in orthogonal coordinates 


Greek Symbols 

3 Crank-Nicolson factor; Roberts’ transformation estimate of 

boundary layer 
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Y 


Ratio of specific heats; angle of coordinate relative to 
the horiztonal 

A Change 

6 Finite difference operator; change; thickness 

K von Karman constant (nominal value 0.43) 

(|) Column vector of dependent variables 

y Viscosity 

p Density 

0 ^ Constant in damping factor (Eq. 11) 

T Shear stress 

ip Mach angle relative to coordinate system 


Subscripts 

b 

lam 

T 


w 

X 

y 

z 

1 

2 

3 


Boundary layer 
Laminar or molecular 
Turbulent; transformation 
Wall 

Associated with the x-direction 
Associated with the y-direction 
Associated with the z-direction 
Associated with the x-direction 
Associated with the y-direction 
Associated with the z-direction 
Free stream conditions 


Superscripts 

n Evaluated at the axial location 

* Value after the first sweep 

** Value after the second sweep 

Ensemble average 
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GOVERNING EQUATIONS 

In this study, calculation of the flow in inlets of rectangular cross 
section is obtained by spatial marching of the equations for conservation of 
mass, momentum and energy. For the class of geometries considered in this 
study (Fig. 1) , the equations are cast in a curvilinear orthogonal coordinate 
system which is aligned with the flow in such a manner that one coordinate 
direction can be identified as the 'primary flow’ direction (Ref. 2). The 
remaining two coordinate directions are used to define the secondary flow. 

The coordinates are chosen as x, y and z with corresponding metric coefficients 
h^, and h^ such that the incremental distance 5s is defined as 

(8s)2 . (h|8x)2 +(hg8y)* + (hjSz)^ a) 


The x~direction is defined as the primary flow or streaim'fise direction while 
y and z define the coordinates in the secondary flow plane at a given stream- 
wise location. The class of geometries considered in this study (rectangular 
inlets) yields metric coefficients which are function of x and y only. 

The governing equations are derived from the Navier-Stokes equations for 
a viscous, heat conducting perfect gas. All variables are ensemble averaged 
and thus the equations apply to both laminar and turbulent flow. Viscous 
terms are simplified by assuming thin boundary layers and all primary flow 
direction (x) diffusion terms are neglected. The equations are as follows: 
Continuity Equation 

-^(h^hjpu) + -^(h,h3^v) + -^(h,h^^w) =0 c 
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Primary Plow (x) Direction Momentum Equation 


r d I ^^^1 1 9 

+ [ IT ‘ ■ ”*ir p'^ 
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( 3 ) 


First Cross Flow (y) Direction Momentum Equation 
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Second Cross Flow (z) Direction Momentum Equation 
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I 1 r d I dh- ■! 

+ [-ay- 

d * dh, , ah, , 

+ -ar > ■ '’2“ai'/“‘' " ^•"aT^'' 
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•ip 
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hjhg dw \ 
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k ^ ay j 

*T — 

3 

az V 

h, ^ arj 


( 5 ) 


It is assumed that the energy equation is represented by the equation 

= constant (6) 


In the above u, v and w represent the mass weighted mean velocity components 
in the x, y and z directions, respectively, e.g., 


u » 

p 


( 7 ) 


where the bar denotes an ensemble average. Also p and p are the ensemble- 
averaged density and static pressure, respectively, p is the effective 
viscosity (defined as the sum of the molecular viscosity and the turbulent 
Boussinesq viscosity, i.e., p = ^lam ^ ^o stagnation enthalpy. 

The pressure, density, velocities and stagnation enthalpy are related to each 
other through the perfect gas equation of state, viz., 

p . 71^ ^[hj,--^(u* + v* +W^)] (8) 
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where y is the ratio of specific heats (assigned to be constant) . The turbulent 
viscosity, is calculated from a mixing length (algebraic) formulation, i.e.. 




m 



1 \* / 

1 du 


1/2 

\ ^ «y j V 

h* ay j [ 

117 IT 

/ 1 h, aw . 



(9) 


where £ is an a priori specified mixing length distribution. In this study 
m — 

the mixing length of McDonald and Camarata (Ref. 3) was utilized. 


i^(y) » 0.09Sjjtanh[ Acy /(0-09Sjj)j^ 


( 10 ) 


where 6, is the local boundary layer thickness, k is the von Karman constant 
b 

(= 0.43), y^ is the distance from the nearest wall, and Ib, the sublayer damping 
factor, is defined by 

^ ■ [ pCy"*"- y'*')] /cr, 

where p(y+ - y+) is the normal probability function, 

y+. 

where t is the local wall shear stress, y"*^= 23 and = 8. 

The partial differential equations (Eqs. (2) through (5)) along with the 
auxiliary relationships (Eqs. (6) and (8) through (10)) represent the set of 
equations which are supposed to describe the flow phenomena inside a rectan- 
gular supersonic inlet. These equations, along with the solution initial values 
at the first streamwise station and the appropriate boundary conditions, will 
be used as the base set of equations which are to be solved in the remainder 
of this report. However, before proceeding further, it is important to 
emphasize that the numerical procedure used in the present work can treat more 
general or alternate forms of the specified auxiliary relationships and 


( 11 ) 


( 12 ) 
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governing equations. For instance the form of the mixing length or the perfect 
gas equation of state could be changed. Additional curvature or stress terms 
(as long as they are not in the primary flow direction) could easily be 
incorporated. 
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BRANCHING SOLUTIONS TO THE GOVERNING EQUATIONS 


At high speed operating conditions the flow upstream of the inlet throat 
is supersonic except for the regions near the walls and thus the governing 
equations in the supersonic core region are essentially spatially hyperbolic 
in nature. In typical supersonic inlets, boundary layer bleed is used to 
minimize distortion and pressure loss generated at incident shock wave 
boundary interactions. The bleed also has the effect of suppressing the extent 
of the subsonic region in the boundary layers, where the differential equations 
are essentially elliptic in nature and downstream influences can effect the 
upstream flow. With sufficient bleed the elliptic (subsonic) regions can be 
kept negligibly small, and if modifications are made to the solution procedure 
in the subsonic region, the governing equations in the supersonic inlet can 
still be solved by a forward marching procedure. 

If the initial value approximations to the Navier-Stokes equations are 
used directly, i.e., without modification, for the case of the rectangular 
supersonic inlet a phenomenon known as ’branching’ generally occurs. 

Basically this phenomenon occurs because of the interaction of the subsonic 
portion of the flow with the supersonic portion of the flow. When an adverse 
pressure gradient occurs (as would happen with the shock boundary layer 
interaction), the subsonic region reacts to this by expanding (see e.g., 
Shapiro, Ref. 4). This in turn causes the supersonic portion of the flow to 
also turn which, of course, causes a further rise in pressure. Since there 
is no damping mechanism for this phenomenon (since streamwlse diffusion terms 
have been removed to create the initial value approximations to the Navier- 
Stokes equations) the result can be an unstable situation. An excellent 
review and interpretation of the branching phenomenon is given by Weinbaum 
and Garvine (Ref. 5). 

Typically the problem of branching with initial value approximations to 
the Navier-Stokes equations has been avoided by explicitly specifying the pres- 
sure, as for Instance is done in classical boundary layer theory. Previous 
work aimed at eliminating this branching problem is largely based on the early 
work of Rudman and Rubin (Ref. 6). Rudman and Rubin solved the equations for 
hypersonic viscous flow over slender bodies with sharp leading edges. In this 
work the streamwlse pressure gradient was neglected and thus the branching 
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phenomenon was suppressed. Based on this two-dimensional work a series of 
three-dimensional problems have been attacked (Refs. 7 to 15). All of these 
three-dimensional solution procedures use iterative techniques and in each case 
the streamwise pressure gradient was neglected. Lubard and Helliwell (Ref. 12) 
also tried setting the streamwise pressure gradient to the explicitly computed 
value but because of the small streamwise pressure gradient associated with the 
geometry, no significant difference was noted over the case with the zero 
pressure gradient assumption. The neglecting of the streamwise pressure 
gradient to suppress branching therefore limited these earlier works to flow 
situations where little streamwise pressure gradient is generated, e.g., flow 
over a cone. Since the purpose of an inlet is to efficiently diffuse the 
external flow for use by the engine by generating a streamwise pressure gradient, 
numerical techniques which neglect the streamwise pressure gradient are not 
sufficient for the analysis of flow in an inlet. The remainder of this section 
will discuss a technique for suppressing branching without neglecting the 
streamwise pressure gradient. 

The technique for implicitly including the effect of the streamwise 
pressure gradient in the initial value approximations to the Navier-Stokes 
equations can be divided Into three parts. First, the equations are assumed 
to reduce to the more conventional boundary layer equations in the subsonic 
regions of the flow. Second, the wall tangency condition is enforced at the 
interface between a supersonic and subsonic region. Finally, the subsonic 
region is kept from directly Influencing the supersonic regions by the use 
of ’type dependent' finite difference operators. The first two parts will 
now be discussed in detail; discussion of the 'type dependent' differencing 
will be deferred until the next section. 

Figure 2 shows the division of a planar cross section (at a given streamwise 
station) into supersonic and subsonic regions for a typical rectangular inlet. 

In this study the initial value approximations to the Navier-Stokes equations, 
represented by Eqs. (2) through (5), were used directly for the supersonic 
region. Since these equations are essentially spatially hyperbolic in behavior 
in the supersonic region, these equations can be replaced by their finite 
difference analogues without further modification and then solved by a forward 
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inarching procedure. In the (thin) subsonic region Eqs. (2) through (5) are 
modified. First it is assumed that in the entire (thin) subsonic region 
the momentum equation normal to the nearest wall is replaced by the conven- 
tional boundary layer approximation that the pressure gradient normal to 
the wall is balanced by the curvature terms. In a region close to a 'y wall' 
the cross flow y-momentum equation is replaced by 


h,h 


3 dy 


ay 


pw 


+ h 


ah 


3 dy 




(13) 


while near a 'z wall' the cross flow z-momentum equation is replaced by 


dp 


ah, 


2 dz 


+ h, 


dh^ 

dz 




(14) 


For the case of rectangular inlets h^ is equal to unity and the metrics are 
functions of x and y only. Hence Eqs. (13) and (14) can be reduced to their 
familiar boundary layer form 


ap I ah 


(15) 


and 



(16) 


The second approximation is termed the wall tangency condition and this 
condition is used to replace the continuity equation on the sonic line by 
the approximation that the flow is parallel to the nearest wall. The tangency 
condition thus becomes a problem of geometry only (Fig. 3). In a region close 
to a 'y wall' (with no bleed) this becomes 

» ton(/^“ X) (17) 
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and in a region close to a ’z wall’ (with no bleed) this becomes (since the 
y’s are a function of y only) 


W = 0 


( 18 ) 


It is to be noted that the normal pressure gradient-curvature relationship 
replaces the cross flow momentum equation everywhere in the subsonic region 
while the wall tangency condition is applied only at the sonic line. In the 
subsonic region the forms of the primary flow or x-momentum equation and the 
continuity equation (below the sonic point) remain unchanged. In addition 
the remaining momentum equation, i.e., the z-momentum equation near a 'y wall’ 
or the y-momentum equation near a ’z wall’, is also unchanged. This last 
feature allows for the existence of y-direction pressure gradients near a 
'z wall’ and z-direction pressure gradients near a 'y wall’. Thus, for 
Instance, glancing shocks on the sidewalls are allowed to generate subsonic 
cross flows as a result of the cross flow pressure gradient generated in the 
supersonic region. The above formulation allows the streamwise pressure 
gradient term to be included without causing the initial value approximate form 
of the Navier-Stokes to exhibit branching. The pressure gradient-curvature 
relationship, along with the wall tangency condition and the appropriate 
differencing, has the effect of allowing the supersonic region to impress its 
pressure on the subsonic region while preventing the subsonic region from 
directly Interacting with the supersonic region. The effect of the subsonic 
region on the supersonic region is indirect and occurs only as a result of 
stream tube expansion, not as the result of rate of change of stream tube 
expansion which it would if the normal velocity were not constrained by the 
wall tangency approximation. This limited interaction between the subsonic 
and supersonic regions is carried over into the numerical difference molecule 
which is ’ type-dependent ’ in the manner in widespread use for transonic 
inviscid calculations (Jameson, Ref. 16). The type-dependent differencing 
of the dependent variables will be discussed in the latter part of the next 
section. 
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NUMERICAL TECHNIQUE 


In this section a brief discussion of the numerical procedure will be 
presented. For more details, the reader is referred to Refs. 2 and 17. 

The general approach is to take the governing system of equations which describe 

the flow in a supersonic inlet (Eqs. (2) through (5) and (15) through (18) 
as previously described) and replace them with finite difference analogues. 

The result of the replacement is a set of coupled nonlinear algebraic equations. 
These equations are linearized by Taylor series expansions about the solution 
at the most recent streamwise station. This results in a series of coupled 
implicit linear equations for the unknown dependent variables at the next 
streamwise station. The terms are then grouped by the coordinate direction 
and then solved by a consistently split, block implicit technique. The 
splitting technique is used to reduce the multidimensional (in this case 
two-dimensional) linearized difference equations to a sequence of one- 
dimensional matrix equations. For three point or less difference formulations, 
these linear one-dimensional difference equations can be written in a block- 
tridiagonal matrix form and solved efficiently and without iteration by 
standard block elimination techniques. To demonstrate how the above formation 
is used, the method will be applied to the continuity equation, Eq. (2). 

The first step is to replace the partial differential equation by its finite 
difference analogue. This results in, for variable centering in the streamwise 
direction, 


+ ( I - [ Sy ( h|h^/3v) + §2 ^ 


n+i 


(19) 


where Ax is the grid spacing in the marching (x) direction, 6^ and 6^ are the 
difference operators in the y and z cross plane directions, respectively, and 
the parameter 3 has been introduced to permit a variable centering of the 
scheme in the x marching direction. If 3 = 1, a fully implicit scheme results 
and 3=1/2 yields the Crank-Nicolson formulation, n refers to the streamwise 
(x) station where the values of the dependent variables are known while 
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n + 1 refers to the next streaTHwise station where the values of the dependent 

variables are unknown. Equation (19) is written at each grid point in the 

plane perpendicular to the streamwise direction, i.e., the y-z plane at x 

station n + 1. Equation (19) is nonlinear in the unknown dependent variables 

and is obviously coupled to the other governing equations through u, v, w and 

til 

p. Linearizing Eq. (19) about the n streamwise station by Taylor series 
expansion yields 

^ n+i^ n+i ^ U + u 2.p U ) n n '^2 ^3 ”2 "s J 

^ P “ 

+ /38y[ h|’*'h;*'(/v'’*' + vV*' - 

+ yS8,[h,''*'h;*' (/«"*' - /w")] 

+ i\- ‘ 0 


Equation (20) is a linear equation in u”"*"^, v^^^, w’^'*'^ and This proce- 

dure can be utilized to derive linear implicit difference approximations to 
the momentum equations (or to the subsonic approximations to these equations) . 
Pressure terms are easily linearized using the equation of state and the perfect 
gas law, for example 



becomes (with constant y and h ) 

o 

(irf ^ +w"w-'] 

+ [ + w^] I 


(21) 


( 22 ) 
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When each of the governing equations is linearized as above, the resulting 
difference approximations can be grouped by coordinate direction and written 
in compact linear matrix difference operator notation as 


A 

Ax 


(<I> 


n+i 


- 


Dy<I> 


n+i 




n+i 


+ S 


(23) 


where <I> is a column vector containing the dependent variables u, v, w and p 
at each grid point in the y-z plane and A is a square 4x4 submatrix. The D 

y 

and are square 4x4 submatrices (written at each grid point in the y-z plane) 

containing elements which are themselves spatial difference operators for the 

y and z directions, respectively. S is a column vector reserved for any 

source terms present. The matrices A, D , D and the column vector S contain 

y ^ n+1 

only quantities which are known, hence Eq. (23) is linear in $ 

Equation (23) could be solved by a direct method but only at considerable 
cost considering the number of grid points required for acceptable accuracy. 

In this study to reduce the computational labor of solving the system of 
equations the Douglas-Gunn splitting technique was utilized. For the case 
considered this results in the equations 

* (<!)* - i") . Dy<t* + Dj*" + S (24) 


for the first sweep and 


-— ( <j) -<!))= + D-O + S 

Ax y 2 


(25) 


** n+1 2 

for the second sweep, where $ is the approximation to $ to order Ax . 

Rather than solve Eq. (25) directly, some computational labor can be saved 

by subracting Eq. (24) from Eq. (25) and solving for the difference across a 

step, i. e. , 


A 

Ax 


( -<I> ) * 


-<P ) 


(26) 
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Additional computational labor can also be saved by observing that during the 
first sweep the implicit coupling between the z-momentum equation and the 
remaining equations is weak and that during the second sweep the implicit 
coupling between the y-momentum equation and the remaining equations is also 
weak. This uncoupling allows one to solve Eqs. (24) and (26) in two parts, 
the first with 3x3 submatrices and the second with 1x1 submatrices rather than 
in one part with 4x4 submatrices. For this system of equations this results 
in approximately 50% savings in computer time in the block elimination proce- 
dure. This concept is termed reducible sub-blocks and is discussed in more 
detail in Ref. 18. 

Before proceeding to a discussion of the forms of the cross plane 
difference operators (6 and 6 ) as applied to the case of the rectangular 

y 2 

inlet, two further areas must be covered, initial conditions and boundary 
conditions. In order to apply Eqs. (24) and (26) it is necessary that all 
values of the dependent variables, viz., u, v, w and p be provided at each 
grid point in the initial y-z plane. Once this is done Eqs. (24) and (26) 
can be marched in the primary flow direction (x) if appropriate boundary 
conditions are specified at the four boundaries corresponding to the 
boundaries of the y-z computational domain. Boundary conditions must be 
specified for each of the four partial differential equations. However 
because of use of the block matrix technique of solving all equations simul- 
taneously (rather than the use of a sequential iterative solution technique) 
coupled boundary conditions can be specified. For instance, one possible 
boundary condition could be to specify that the pressure normal to the cowl 
is zero , i . e . , 


dy 


= 0 


(27) 


However the pressure has been eliminated from the system by use of the perfect 
gas law and the assumption of constant stagnation enthalpy. Thus Eq. (27) 
becomes after linearization, a coupled relationship, Eq. (22) between the 
dependent variables. However the virtue of the block implicit approach is 
that this coupled but linearized relationship is easily implemented as a 
boundary condition. 
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Turning to the problem of 'type dependent' differencing, the form of the 
difference operators utilized In the present study Is a necessary part of 
the technique of prohibiting the subsonic flow from causing the solution to 
branch. In the supersonic portion of the flow all cross plane derivatives 
are represented by means of standard three point central difference formulae. 

In the subsonic portion of the flow all pressure derivatives in the direction 
normal to the nearest wall are taken as two point one-sided away from the wall. 
This allows the subsonic pressure to be implicitly set by the supersonic flow. 
In addition, in the subsonic region the convective derivatives normal to the 
nearest wall are differenced as two point one-sided towards the wall. Finally 
on the sonic line, the convective derivatives normal to the nearest wall are 
differenced as two point one-sided away from the wall. 


20 



TEST CASES 


In order to evaluate the computational procedure described in the 
previous sections, two test cases were run with the PEPSIS computer code. 

The first case considered was that of a Mach number 2.5 flow entering a 
square duct. At some point downstream of the entrance air was blown into 
the duct through both the bottom and top walls. Further downstream the 
blowing was terminated. The intent of this test case was to determine if 
the PEPSIS computer code would correctly predict the expansion and compression 
waves produced by the initiation and termination of blowing. The second 
case considered was the flow in a Mach number 3.0 high-performance, ’’two- 
dimensional" mixed-compression, inlet system (Refs. 19 to 21). Although 
termed "two-dimensional" the inlet had a rectangular cross-section and the 
computations reflected this three-dimensionality. For the inlet case 
experimental data were available and hence comparisons were made between the 
predicted and measured results. Details of the two test cases run, the 
results of the theoretical predictions, and (for the inlet case) comparisons 
between theory and measurements will next be discussed in aetail. 

The first test considered was a uniform flow entering a square duct 
(Figure 4). The free stream Mach number was 2.5, while the Reynolds number 
based on free stream properties and the duct height was 5000. The viscosity 
was asstimed to be constant. The primary direction of the flow was parallel 
to the walls of the duct. From 0.36 to 0.44 duct heights downstream of the 
duct entrance, air was normally injected into the oncoming flow from both 
the floor and top surface of the duct. The injection rate was varied in the 
spanwise plane as is shown in Figure 5 . For this case a 10 x 10 uniform 
finite difference grid was utilized. Because of the uniform nature of the 
flow entering the square duct and the geometry of the duct two planes of 
symmetry can be associated with this problem, and hence the flow need only 
be solved in the quarter plane. On the planes of symmetry, symmetry boundary 
conditions were utilized, while on the walls the three no— slip conditions were 
used for the velocity components and the normal momentum equation was used as 
the fourth condition. Figures 6 through 8 show the transverse pressure 
coefficient distributions (at four spanwise stations each) at three selected 
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streamwise stations, l.e., at streamwlse locations 0.522, 0.810, and 2.213 
duct heights from the duct entrance. 

In Figure 5, which is slightly downstream of the cessation of blowing, 
the compression disturbance can be seen propagating across the duct. Since 
the amount of blowing decreases as the sidewall Is approached, the strength 
of the disturbance is smaller in the region of the sidewall. In Fig. 6 the 
effect of the expansion (caused by the cessation of blowing) can be seen to 
have somewhat cancelled the effect of the compression in the region of the 
wall, and by the time the flow has reached the last streamwise station (Fig. 6) 
the interaction effects of the compression and expansion have to a great 
degree cancelled each other out. Figure 9 shows the streamwise variation of 
the pressure coefficient at selected transverse locations on the vertical 
plane of symmetry (curves (1) , (2) and (3) correspond to transverse locations 
of 0.05, 0.22 and 0.45 duct heights respectively). From careful examination 
of this plot it can be deduced that over the streamwise domain of this calcu- 
lation the compression wave has reflected off the horizontal plane of symmetry 
and has propagated back down to the bottom wall. On the other hand, the 
expansion wave has propagated upward, reflected off the plane of symmetry, 
but not yet reached the bottom wall. 

For the second test case, the PEPSIS computer code was used to calculate 
the flow field for the NASA Ames Research Center large-scale variable-geometry 
inlet system operating at a design Mach number of 3.0 with a Reynolds number 
per foot of 2.2 x 10^. A schematic of the inlet system is shown in Fig. 10. 

For more details of this experimental apparatus the reader is referred to 
Refs. 19 to 21. The inlet capture area is 14 inches square with the ramp 
initially inclined 7° to the horizontal. A variable-ramp drive mechanism 
allows the second portion of the ramp to pivot around a point 28 inches 
from the ramp leading edge. The design condition (and that used for the 
computations) was a second ramp angle inclined at 14° to the horizontal while 
the cowl surface is parallel to the horizontal. The coordinate system used 
for the computation of the flow field is shown in Fig. 11. This coordinate 
system was generated by Charles Towne of NASA Lewis Research Center with the 
ADD computer code (Ref. 22). The ADD code uses a Schwartz-Cristof f el trans- 
formation technique to calculate the conformal coordinate system (hj^ = h^) 
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associated with a given geometry. The ADD code will not generate a coordinate 
system that allows for a step discontinuity in curvature (as happens when the 
ramp angle discontinuously changes from 7° to 14°, but instead generates a 
smooth curve rather than a discontinuity. For instance, inspection of Fig. 11 
shown that the initial ramp angle is not 7°, but instead transitions from 
0° to 7° gradually. The second ramp angle, 14°, also does not occur dis- 
continuously . 

The coordinate system shown in Fig. 11 was used by the PEPSIS computer 
code to calculate the flow field inside the NASA Ames inlet . It was assumed 
that the initial flow was at uniform free stream conditions. For the first 
runs on the NASA Lewis Research Center's UNIVAC 1100 a 20 x 20 nonuniform 
mesh was used in the cross-sectional plane with (as can be seen in Fig. 11) 
packing of grid points in the region of the ramp and cowl. The PEPSIS 
computer code was marched 181 streamwise steps, the computation terminating 
just upstream of the inlet throat. For the first 10 streamwise stations a 
marching step size. Ax, of 0.025 was used. From streamwise station 10 to 20, 
the marching step size was increased from 0.025 to 0.04. After station 20 
the step size remained constant at 0.04. The cowl was encountered at stream- 
wise station 70. (It is to be noted that the streamwise grid of Fig. 11 only 
plots every other marching station.) Ax is the computational step size and 
not the physical step size, which is given by h^(x)Ax, thus accounting for 
the varying physical step size when Ax is a constant. 

The boundary conditions utilized for the computation in the inlet can be 
divided into three categories: (1) boundary conditions on walls, (2) boundary 

conditions on the computational surface upstream of cowl, and (3) boundary 
conditions on the plane of symmetry. The wall boundary conditions require 
that the two components of the velocity parallel to a wall obey the wall 
function boundary condition (see Appendix B) and that the velocity component 
normal to the wall be zero (or equal to the bleed rate) . In addition the 
momentum equation in the direction normal to a wall is presumed to be satis- 
fied. On the computational surface upstream of the cowl Mach line extrapol- 
ation was utilized (see Appendix B) . Initially some difficulty was encountered 
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with these boundary conditions. Preliminary analysis indicates that this 
problem was associated mainly with the poor accuracy of the metric infor- 
mation in the vicinity of the boundaries. The last category of boundary 
conditions was applied on the transverse plane of symmetry. This computa- 
tional plane of symmetry exists only because of the geometry of the inlet 
considered and because a 0° yaw case was considered. The symmetry conditions 
require that the transverse derivatives of the streamwise velocity, u, 
transverse velocity, v, and the density, p, be zero at the plane of symmetry. 
In addition, the spanwise velocity, w, must be zero. Because of the magni- 
tude of the Reynolds number of the flow inside the inlet, the entire flow 
was presumed to be turbulent. For the case considered the previously 
discussed mixing length model of McDonald and Camarata (Ref. 3) was utilized. 
Finally, for this calculation the effect of bleed was neglected. 

Results of the computation are presented in Figs. 12 to 17. Figure 12 
shows the inviscid shock pattern. The pattern for the above case was com- 
puted by two-dimensional method of characteristics (Ref. 19). Figure 13 
shows the "shock” pattern as calculated by the PEPSIS computer code. Except 
for the shock produced by the discontinuous 7° -14° bend in the ramp all of 
the inviscid shock waves are reproduced by the PEPSIS calculation. It is 
likely that the absence of the shock due to the ramp discontinuity is caused 
by the previously discussed smoothing produced by the ADD code coordinate 
system. This smoothing causes the PEPSIS code to produce a series of 
compression waves rather than an easily identifiable shock wave. Figures 
14 and 15 show the streamwise distribution of pressure along both the ramp 
and cowl at the plane of symmetry. It is evident that the PEPSIS computations 
agree very well with the experimental results, both in the positioning of 
and the strength of the shock waves. Figures 16 and 17 compare the two- 
dimensional method of characteristics ramp and cowl pressure distributions 
with those produced by the PEPSIS computer code. The PEPSIS results compare 
favorably with the method of characteristics solution. The only difference 
is that the PEPSIS results are slightly more diffuse than the method of 
characteristics solution. 
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USERS' MANUAL 


The PEPSIS users' manual is meant to serve as a guide in helping the 
user perform successful runs with the PEPSIS computer program. The manual 
is divided into six parts. The first part is a flow diagram to help the 
user understand the basic flow of the program. Because of the size of the 
code, a detailed flow diagram would be rather cumbersome to use and would 
probably by of little use to the user. Therefore the flow diagram provides 
only the 'big picture' of the code. If a user is interested in the details, 
it is suggested that he consult the program listing. The second part describes 
in detail the input required to run a case with the PEPSIS computer code. 

The vast majority of the input is input via the NAMELIST format. There are 
two advantages to the NAMELIST format, viz., that if the default values are 
acceptable, the user need not input that variable and that, within a given 
NAMELIST, the order in which the variables are entered is irrelevant. The 
third part of the users' manual describes the error conditions which can occur 
in PEPSIS and the corrective action which should be taken if they occur. 

The fourth part describes the PEPSIS FORTRAN variables. The FORTRAN variable 
(with dimensions for dimensional variables) , their COMMON block for the vari- 
able, and the description of the variable is presented for the variables in 
alphabetical order. The fifth part describes the storage devices which are 
required to make a run with PEPSIS. If general orthogonal coordinates are 
used, the present version of PEPSIS reads data off unit 10 by means of binary 
reads. The main reason that this was done was to economize on computer run 
time as formatted reads take considerably more time than binary reads. 

However if the metric data generated on unit 10 is done on say a UNIVAC com- 
puter and the PEPSIS run is to be made on say a CDC computer, the binary 
reads must be changed to formatted reads because of the incompatability 
between UNIVAC binary and CDC binary. The required changes must be made in 
SUBROUTINE GEOTRB. Finally, the sixth part of the users' manual gives the 
input and output for a typical PEPSIS run. 

Before proceeding to the first part, it is perhaps wise to give some 
general background that will help in interpreting that which follows. 

Figure 18 is a schematic of a typical cross flow section at a given axial 
station. In the PEPSIS computer code, the transverse direction is denoted 
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as the y-direction while the spanwise direction is denoted as the z-direction. 
The lower boundary is called the 1 surface, the upper boundary the 2 surface, 
the left boundary the 3 surface and the right boundary is called the 4 
surface. The PEPSIS logic assumes that the ramp is on a 1 surface and that 
the cowl is on a 2 surface. The 3 surface is assumed to be a plane of 
symmetry while the sidewall must be on a 4 surface. 
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PEP SIS Input 


Except for a title card and two input variables which are needed by the 
plot file, all PEPSIS input is input through NAMELISTS $LIST2 through $LIST8. 
The input variables in each of the NAMELISTS will be described and if a default 
value exists, its value will be given. 

Plot File Input 


Card 1 

Columns 

Format 

Variable 

Function 

1-32 

5A6,A2 

TITLE(I) 

Title Card 

Card 2 

Columns 

Format 

Variable 

Function 

1-2 

12 

ISYM 

Reciprical of Symmetry 

3-12 

FIO.O 

SYSTEM 

SYSTEM=1 - Quasi-Cartesian 
Coordinates 

SYSTEM=2 - Quasi-Cylindrical 
Coordinates 


Namelist Input Description 


Namelist or 
variable name 

$LIST2 

UZERO 

YZERO 

TZERO 

VISCOS 

RG 

USCALE 

IWR 


Description 
Reference Variables 

Reference value of velocity (ft/sec). 

Reference length (ft) . 

Reference temperature (°R) , Default 416.0 °R. 

2 

Reference kinematic viscosity (ft /sec). 

Gas constant (ft-lbf /slug-°R) , Default 1716.3 
Scale factor for metric information from ADD code, 
Flag to print out namelists. Default value is 1. 
IWR=1. Print namelists. 

IWR=0. Do not print namelists. 
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IWR=0. 



Namelist or 
variable name 


Description 


$LIST3 

lAXI 


ICAR 


$LIST4 

YS(1,1) 

YS(2,1) 

YSCl,2) 

YS(2,2) 

NE(IADI) 

IGRID(IADI) 


Geometry Options 

Axisymmetric geometry sentinel. Default value is 0. 

lAXI = 0. No axisyrmnetric option. 

lAXI = 1. Axisymmetric option used. 

Cartesian geometry sentinel. Default value is 1. 

ICAR =0. No Cartesian option. 

ICAR = 1. Cartesian geometry used. 

If ICAR = 0 and lAXI =0, general orthogonal 
coordinates are used. For this option storage unit 
10 must be generated from the ADD code (or equivalent) 
(Ref. 22) . 

Geometry and Restart Options 

Nondimensional value of y at the number 1 surface. 

Nondimens ional value of y at the number 2 surface. 

Nondimensional value of z at the number 3 surface. 

Nondimensional value of z at the number 4 surface. 

Number of grid points in the y(IADI = 1) and 
z(IADI = 2) directions. 

Roberts’ grid transformation sentinel for y and z 
directions, respectively. Default values are 1. 

(See Appendix A). 

IGRID(IADI) = 0. Equally spaced grid points. 

IGRID(IADI) = 1. One boundary layer grid at surface 
number 1 (lADI = 1) or 3 (lADI = 2) . 

IGRID(IADI) = 2. Boundary layers at both surface 
numbers 1 and 2 (lADI = 1) or surface numbers 3 and 4 
(lADI = 2). 

IGRID(IADI) = 3. One boundary layer grid at surface 
number 2 (lADI = 1) or 4 (lADI = 2). 
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IIIIIH I ■III! II 


Namelist or 
variable name 

EPS(IADI) 


NS 

X1,X2,XENTR,AP 


IRSTIN 


IRSTOT 

JRSTIN 

JRSTOT 

NFILE 


Description 

Estimate of boundary layer thickness for use with 
Roberts' grid transformation. If IGRID(IADI) = 2, 
the value of EPS(IADI) Is the same on both surfaces. 
EPS(IADI) equals the percentage of the distance 
YS(2,IADI) - YS(1,IADI). EPS corresponds to the 3 
used In Appendix A. 

Number of stations to be marched during a run. 

X2-X1 is the nominal value of AX in the equation 

X(I) = X(I-l) + AP*(X(I-1) - X(I-2)) 

where X(ISTART) = XENTR 

and XCISTART+I) = X(ISTART) + AX 

and I varies from ISTART+2 to NS. For a restart 

ISTART=IRSTIN while for a new case ISTART=1, AP 

represents the rate at which the marching step size 

grows (AP > 1) or diminishes (AP < 1) . Default values 

of XENTR and AP are 0.0 and 1.0, respectively. The 

value of XENTR must be input during each restart 

run and must equal the value of X at the restart 

station. 

Marching station number where data is to be read 
in for restart case. Default value is 0. 

IRSTIN = 0. Not a restart case. 

IRSTIN ^ 0. Restart case started from data at 
station IRSTIN. 

Interval for saving restart information. 

Default value is 0. 

IRSTOT =0. No restart information is saved. 

IRSTOT 4 ^ 0. Information is saved at each IRSTOT^h 
station. 

Logical file name from which restart information is to 
be read. Default value is 11. 

Logical file name onto which restart information is to 
be written. Default value is 11. 

Sequence number in JRSTIN of the desired restart 
information. 
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Namelist or 
variable name 

NSAVED 

AMDOTD(I,JX,N) 

XJXTW 

$ LIST5 

ILA14 

^IJ_ST6 

YSLOT(IADI) 

NBC0N(IVAR,N,1) 


Description 

Number of restart stations saved on JRSTOT . Default 
value is 0. On a restart by setting JRSTOT = JRSTIN 
and NFILE = NSAVED one file can be used for both read 
and write without destroying the information 
previously saved. 

Distribution of nondimensional bleed velocity at 
marching station JX and surface N. The index I 
represents the spanwise (for N = 1 and 2) or transverse 
(for N = 3 and 4) grid locations for the distributions. 
Default values are 0.0 AMDOTD(20,200,4) . 

Marching direction location at which the cowl begins. 
Default value is 0.0. 

Turbulent Option 

Sentinel for laminar flow. Default value is 1. 

ILAM = 0. Turbulent flow. 

ILAM = 1. Laminar flow. 

Mixing Length Constant 

Representative location of outer edge of boundary 
layer, 6^ of Eq. (10). 


Boundary Conditions 

Boundary condition of the IVAR^^ variable on the 
surface. 

X-Momentum Equation (NVl = 1) 

NBC0N(NV1,N,1) =0. U set equal to 0.0. 

NBC0N(NV1,N,1) = 1. First derivative of u set equal 
to 0.0 when N ^ 2. For N = 2, u is set by Mach line 
extrapolation in region upstream of the cowl. (See 
Appendix B) . Downstream of the cowl appropriate wall 
boundary condition is set automatically by code. 

NBC0N(NV1,N,1) = 2. Second derivative of u set equal 
to zero. 

NBC0N(NV1,N,1) = 3. Not used. 
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Namelist or 
variable name 


Description 

NBC0N(NV1,N,1) = 4. Value of u set by wall functions. 
See write-up on NBC0N(NV2,N,1) = 4 and NBC0N(NV3,N,1) 

= 4. (See Appendix B) . 

Y-Momentum Equation (NV2 = 2) 

NBC0N(NV2,N,1) = 0. V set equal to AMD0TD(IZ, JX,N) 
where IZ is the spanwise grid location, JX is the 
marching location and N is 1 or 2. For N equal to 
3 or 4, V set equal to 0.0. 


NBC0N(NV2,N,1) = 1. First derivative of V set equal 
to 0.0 when N?^2. IfN = 2,Vis set by Mach line 
extrapolation upstream of the cowl. Downstream of the 
cowl appropriate wall boundary condition set 
automatically by code. 


NBC0N(NV2,N,1) = 2. Second derivative of V set equal 
to zero. 

NBC0N(NV2,N,1) = 3. Not used. 

NBC0N(NV2,N,1) = 4. Value of V set by wall functions. 
Cannot be used for N = 1 or 2. If used for N = 3 or 4, 
values of NBC0N(NV1,N,1) must be set equal to 4. 


Z-Momentum Equation (NV3 = 4) 


NBC0N(NV3,N,1) =0. W set equal to AMDT0D(IY, JX,N) 
where lY is the transverse grid location, JX is the 
marching location and N is 3 or 4. For N equal to 
1 or 2, W is set equal to 0.0. 


NBC0N(NV3,N,1) = 1. First derivative of W set equal 
to 0.0. 

NBC0N(NV3 ,N, 1) = 2. Second derivative of W set equal 
to 0.0. 


NBC0N(NV3,N,1) = 3. Not used. 

NBC0N(NV3,N,1) = 4. Value of W set by wall functions. 
Cannot be used for N = 3 or 4. If used for N = 1 or 
2, Values of NBC0N(NV1,N,1) must be set equal to 4. 


34 



Namelist or 
variable name 


$LIST8 

IPLOT 



NBCON(NRHO,N, 1) = 0. Pressure specified to be that 
of pressure at previous station. 


NBC0N(NEH0,N,1) = 1. Apply normal momentum equation 
boundary condition when N?^2. IfN=2, pis set 
by Mach line extrapolation upstream of the cowl. 
Downstream of the cowl appropriate wall boundary 
condition set automatically by the code. 


NBC0N(NEH0,N, 1) = 2. Second derivative of density 
set equal to 0.0. 


Plot Information 


Plot option sentinel. Default value is 0. 
IPLOT =0. No plots. 

IPLOT = 1. Plot file created. 
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Error Conditions in PEPS IS 


SUBROUTINE GEOTRB 


In SUBROUTINE GEOTRB there are three modes of failure when using the 
general orthogonal geometry capability. All failure modes indicate that an 
attempt is being made to calculate metric information in a region where the 
ADD computer code (or its equivalent) did not generate data. The error 
message is printed in the dayfile. 


'STOP GEOTRB' 


'STOP GEOTRB 1’ 


'STOP GEOTRB 2' 


SUBROUTINE RESTAR 


Indicates that an attempt is being made to obtain 
metric information at a value of x less than that 
at the first x-station for which data was generated 
by the ADD code. 

Indicates that an attempt is being made to obtain 
metric information at a value of y less than that of 
the first y-station for which data was generated by 
the ADD code. 

Indicates that an attempt is being made to obtain 
metric information at a value of y greater than that 
at the last y-station for which data was generated 
by the ADD code. 


There are several modes by which SUBROUTINE RESTAR can fail. All involve 
improper use of the restart file. A message, 'STOP RESTAR', will be printed 
in the dayfile and a message 'RESTART REQUESTED AT STATION (Station Number) 

BUT STORED INFORMATION AT SEQUENCE (File Number) IS AT STATION' (Station 
Number) . This message appears because the value of JX (marching station number) 
does not match the value of IRSTIN (see input) . This can be caused by either 
an incorrect value of NFILE (see input) and a proper value of IRSTIN or an 
Incorrect value of IRSTIN and a proper value of NFILE or incorrect values of 
both variables. Another failure mode would be for the file JRSTIN to be null 
or for the value of IRSTIN to be incorrect. These failure modes will cause 
the program to cease while in the NTRAN routine. 


SUBROUTINE MINVRS 

The mode of failure in this routine indicates that a matrix is singular 
and thus cannot be inverted. A message, 'STOP MINVRS', is printed in the day- 
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file and a message 'MATRIX SINGULARITY' followed by four integers is printed 
out. This usually indicates that an incorrect boundary condition(s) has been 
chosen. It is suggested that the values of NECON be re-evaluated. 

FAILURE IN SQR T 

If an attempt is made to take the square root of a negative number, a 
failure will occur. Usually this is due to a breakdown of the solution 
procedure when the flow wants to recirculate - (an impossibility with a para- 
mobilized method) . The case to be run should be re-evaluated to determine if 
more bleed should be added to prevent the recirculation. 
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PEPSIS FORTRAN Variables 


FORTRAN 

SYMBOL 

COMMON 

BLOCK 

DESCRIPTION 

A(5f 16«?0) 


STORAGE FOR LINEARIZATION BLOCK MATRIX 
COEFFICIENTS - EOUIVALENCEO TO BLK 

AB 

INPUT 

GAS CONSTANT / SCECIFIC HEAT 

ABEN6 

INPUT 

A6ENG s (GAMHA-1) • CMACH 

A6(20«9*2) 

WEIGHT 

DIFFERENCE WEIGHTS 

AG7 

WEIGHT 

INVERSE OF STREAMWISE STEP SIZE 

AH00T0(20f 200t4) 

CET 

BLEED RATES 

ATS(19*20) 

STORE 

STORAGE ARRAY FOR LINEARIZATION COEFFICIENTS 

B 

INPUT 

B = -0.5 * AB 

BLK<2500) 

BLKMM 

GENERAL PURPOSE STORAGE ARRAY 

BLKl (2500) 

BLKMM 

GENERAL PURPOSE STORAGE ARRAY 

BWD 

BWOIF 

CRANK-NICHOLSON FACTOR 

BWDPl 

BWDIF 

BWDPl = BWD ♦ 1 

CMACH 

INPUT 

MACH NUMBER 

C0ORD(20*4) 

COORD 

GEOMETRIC INFORMATION NEEDED FOR PLOT FILE 

OW 

BCONW 

SOURCE TERM - USED BY BOUNDARY CONDITION 
SUBROUTINE 

OWT (5*4) 

BCONW 

STORAGE FOR LINEARIZATION COEFFICIENTS IN 
BOUNDARY CONDITION SUBROUTINE 

EPS(2) 

INPUTA 

MESH DISTRIBUTION FACTOR 

FF(2.20.3) 

DFNCTN 

GENERAL STORAGE ARRAY - USED MAINLY FOR 
GEOMETRY 

F6(fl*2t20) 

GEOMV 

STORAGE FOR METRIC COEFFICIENTS AND 
DERIVATIVES 

F(l4f 1) 


STORAGE FOR PLOT INFORMATION - EOUIVALENCEO 
TO BLK 

GAMM (20,2) 

GANG 

ANGLE OF COORDINATE LINES RELATIVE TO 
THE HORIZONTAL 

GAMMA 

INPUTA 

RATIO OF SPECIFIC HEATS 

GAMMR 

INPUTA 

PRODUCT OF GAMMA TIMES GAS CONSTANT 

lAOI 

NTRAN2 

ADI SWEEP DIRECTION 

lADSHF (2) 

GENDIM 

SHIFT INDEX 

IJ 

INDEX 

LINES ON WHICH ADI IS BEING PERFORMED 

IJM 

NTRANZ 

COUNTER FOR NTRANNING SUBROUTINES 

lAXI 

PCOR 

SENTINEL FOR AXISYMMTRl’c COORDINATES 

ICAR 

PCOR 

SENTINEL FOR CARTESIAN COORDINATES 

ILAM 

INPUT 

SENTINEL FOR LAMINAR AND TURBULENT FLOW 

10 

RESTRT 

SENTINEL FOR READING AND WRITING 
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PEPSIS FORTRAN Variables 


FORTRAN 

SYMBOL 

COMMON 

BLOCK 

DESCRIPTION 

IPLOT 

INPUTA 

SENTINEL FOR PLOT OPTION 

IRSTIN 

RESTRT 

AXIAL STATION NUMBER FOR RESTART 

IRSTOT 

RESTRT 

INTERVAL FOR RESTART INFORMATION 

ISSC20*2) 

CISS 

LOCATION OF SONIC LINE 

ISSD(20«2) 

CISS 

LOCATION OF SONIC LINE 

ISYM 

BPLOT 

RECIPORICAL OF SYMMETRY 

IWALLF 

BCONW 

SENTINEL which DETERMINES IF WALL FUNCTION 
LOGIC IS NEEDED 

JRSTIN 

RESTRT 

LOGICAL FILE FROM WHICH RESTART INFORMATION 
IS READ 

JRSTOT 

RESTRT 

LOGICAL FILE ON WHICH RESTART INFORMATION IS 
WRITTEN 

JX 

INDEX 

marching STATION INDEX 

K9 

NTRANZ 

STATUS CHECK 

LGP(ItJ) 


LGP(I.J) = 1 ♦ (J-1) * NPOINT - statement 
FUNCTION DEFINING GRID POINT ADDRESS 

LZ 

PARAMl 

LZ = MZVAR * MLEVEL » NPOINT « NZ 

MEZ 

NTRANZ 

MEZ = NOUTER * 1 

MLEVEL 

PARAMl 

MAXIMIUM levels of STORAGE 

MU 

VARNUM 

INDEX FOR LAMINAR VISCOSITY 

MUT 

VARNUM 

INDEX FOR TURBULENT VISCOSITY 

MZ 

PARAMl 

MZ = LZ ♦ 1 

MZVAR 

PARAMl 

MAXIMIUN NUMBER OF STORAGE VARIABLES 

N<7) 

BCONW 

INDICIES USED BY BOUNDARY CONDITION ROUTINE 

NBCON <5f 4» ?) 

BCONW 

INDEX WHICH SPECIFIES TYPE OF BOUNDARY 
CONDITION 

NCAP 

BCONW 

FLAG WHICH TELLS WHEN END CAP LOGIC IS TO 
BE USED 

NFILE 

RESTRT 

SEQUENCE NUMBER OF RESTART INFORMATION 

NCOL 

6ENDIM 

NCOL = 3 * NROW ♦ 1 

NCOLS (2) 

GENDIM 

NCOL (I) = 3 * NROWSm ♦ 1 

NEN6 

VARNUM 

INDEX FOR ENTHALPY 

NEZ 

NTRANZ 

NEZ = NUPPER ♦ 1 

NHl 

VARNUM 

INDEX FOR METRIC COEFFICIENT IN X - DIRECTION 

NH12 

VARNUM 

INDEX FOP DERIVATIVE OF X METRIC IN Y DIRECTION 

NH2 

VARNUM 

INDEX FOR METRIC COEFFICIENT IN Y - DIRECTION 

NH21 

VARNUM 

INDEX FOR derivative OF Y METRIC IN X DIRECTION 
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PEPSIS FORTRAN Variables 


FORTRAN 

COMMON 

DESCRIPTION 

SYMBOL 

BLOCK 


NH3 

VARNUM 

INDEX FOR METRIC COEFFICIENT IN Z - DIRECTION 

NH31 

VARNUM 

INDEX FOR DERIVATIVE OF 7 METRIC IN X DIRECTION 

NH32 

VARNUM 

INDEX FOR DERIVATIVE OF Z METRIC IN Y DIRECTION 

NIM 

NTRANZ 

COUNTER FOR NTRANNIN6 SUBROUTINES 

NIN 

PARAMl 

MAXIMIUM LINES OF STORAGE IN CORE AT ONE TIME 

NLEN 

VARNUM 

INDEX FOR MIXING LENGTH 

NE(2> 

GRID 

NUMBER OF GRID POINTS IN Y AND Z DIRECTIONS 

NOUTER 

NTRANZ 

NOUTER = NE (3-1 ADD - 1 

NPOINT 

PARAMl 

MAXIMIUM NUMBER OF GRID POINTS IN Y OR 
Z DIRECTION 

NREAO 

NTRANZ 

COUNTER FOR NTRANNING SUBROUTINES 

NRHO 

VARNUM 

INDEX FOR DENSITY 

NROW 

GENDIM 

NUMBER OF ROWS IN COEFFICIENT MATRIX 

NROWS(2) 

GENOIM 

NUMBER OF EQUATIONS SOLVED SIMULTANEOUSLY 
IN A SET 

N5 

GRID 

NUMBER OF STATIONS TO BE MARCHED 

NSAVEO 

RESTRT 

NUMBER OF RESTART STATIONS SAVED 

NSETA 

GENDIM 

SENTINEL FOR INDICATING WHICH SET OF 
EQUATIONS BEING SOLVED 

NSLAB 

NTRANZ 

NSLAB s! MLEVEL * MZVAR * NPOINT 

NTAU 

VARNUM 

INDEX FOR STRAIN TENSOR CONTRACTION 

NUPPER 

OUTA 

NUPPER = NE(IADI) - 1 

NVl 

VARNUM 

INDEX FOR X - DIRECTION VELOCITY 

NV2 

VARNUM 

INDEX FOR Y - DIRECTION VELOCITY 

NV3 

VARNUM 

INDEX FOR Z - DIRECTION VELOCITY 

NZ 

PARAMl 

MAXIMIUM NUMBER OF GRID POINTS IN Z 
DIRECTION 

PZERO 

8PL0T 

REFERENCE PRESSURE 

R6 

INPUTA 

GAS CONSTANT 

RMl 

INPUT 

INVERSE OF THE REYNOLDS NUMBER 

RZERO 

INPUT 

REFERENCE DENSITY 

SYSTEM 

BPLOT 

PLOT OPTION 

TITLE<6) 

BPLOT 

TITLE INFORMATION 

TZERO 

INPUT 

REFERENCE TEMPERATURE 

USTAR (20«4) 

USTARR 

FRICTION VELOCITY 

OZERO 

INPUT 

REFERENCE VELOCITY 
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PEPS IS FORTRAN Variables 


FORTRAN 

COMMON 

DESCRIPTION 

SYMBOL 

BLOCK 


VKC 

TURA 

VON KARMAN CONSTANT 

VISCOS 

INPUTA 

LAMINAR KINEMATIC VISCOSITY 

X(S0) 

GRID 

MARCHING STATION LOCATIONS 

XENTR 

INPUT 

MARCHING STATION LOCATION AT BEGINNING OF A 
CALCULATION 

XJXTtf 

CWALL 

AXIAL LOCATION OF COWL 

Y(20t2) 

GRID 

PHYSICAL DISTANCE IN COMPUTATIONAL DIRECTION 

YS<2.2) 

GRID 

NONDIMENSIONAL EXTENTS OF COMPUTATIONAL 
DOMAIN 

YSLOT(2> 

IPROFF 

BOUNDARY LAYER THICKNESSES 

YSAVE{20*2) 

GRID 

COMPUTATIONAL COORDINATES 

YZERO 

INPUT 

REFERENCE LENGTH 

Z(MtNtI*J) 


EQUIVALENT OF ZZ - STATEMENT FUNCTION 

ZBLK(2500) 

BLKM 

GENERAL PURPOSE STORAGE ARRAY 

ZZ ( 17*3*B0) 

VARIAB 

GENERAL PURPOSE STORAGE FOR OEPENDENT 
VARIABLES 

ZZZ(MfN«LGP(I f J) } 


EQUIVALENT OF ZZ - STATEMENT FUNCTION 
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Storage Devices Needed to Run PEPSIS 

The following storage devices are needed to successfully run the PEPSIS 
computer program. 

Unit Number Use 

8 Plot tape, needed only when input value of IPL0T=1. 

This file must be saved after a run. 

9 General purpose NTRAN mass storage device. Assumes 

word addressable device. No need to save this file. 

10 Storage device for metric data generated by the ADD 

computer code. Only needed when IAXI=0 and ICAR=0, 
i.e., for general orthogonal coordinates. 

JRSTIN Logical file from which restart information is to be 

read. Default value is 11. 

JRSTOT Logical file onto which restart information is to be 

written. Default value is 11. This file must be 
saved after a run. 

Input and Results for a Sample Case 

Table I presents the input for a sample case the results of which are 
presented in the following pages. This run is taken from one of the first 
runs made on the inlet test case previously reported. The first part of 
the input is associated with the title and options for the plot option. Input 
in $LIST2 requires that no NAMELIST information will be printed (IWR=0) and 
inputs the reference values that yield a Mach number equals 3.0 flow. The 
$LIST3 data yields the general orthogonal geometry option and thus requires 
that the ADD computer geometry be generated on unit 10. $LIST4 data sets up 
the computational domain, the bleed schedule and restart information. This 
case was run three streamwise stations (to keep the printout manageable 
for the report) . The value of IRSTOT requires that restart information be 
stored each three stations. A value of IRSTIN=0 indicates that this is the 
first run of the computation, i.e., not a restart. $LIST5 datum tells that 
this case is for a turbulent run. The input for $LIST6 contains the estimate 
of the boundary layer thickness while the $LIST7 contains the boundary 
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conditions. These boundary conditions are typical of those required for 
turbulent inlet calculations. It is to be noted that the 2 surface boundary 
conditions are not those required for a cowl. Rather they are the boundary 
conditions associated with Mach line extrapolation. The code will automatically 
generate the cowl boundary conditions when the marching direction distance 
exceeds XJXTW. The final NAMELIST, $LIST8, is void because the default value 
of IPLOT is zero and this is the value desired. It is to be noted that even 
though no values are input in $LIST8 the cards $LIST8 and $END are required. 

The output from PEPSIS for the above input is presented in Table II. 

The first two pages consist of the central difference operators for the 
cross flow directions. Next the metric information at the starting station 
is displayed. The BETA on the far right is the streamwise coordinate angle. 

Next the initial conditions are printed. For turbulent flow this consists of 
the cross flow plane distribution of x-velocity, y-veloclty, z-velocity, density, 
stagnation enthalpy, turbulent viscosity, turbulence kinetic energy, dissipation, 
mixing length, shear, pressure coefficient divided by 2, Mach number and finally 
temperature. All of these variables are in nondimensional form, e.g., veloci- 
ties are nondimensionalized by UZERO. The IZ across the top of the distri- 
butions refers to z grid point location and z refers to the actual z location. 

The integers on the left-hand side refer to the transverse grid point location 
while the numbers to their immediate right are the corresponding values of 
the physical distances (corresponding to the FORTRAN variable Y) and are 
actually the integral of the computational y times h 2 . Above each distri- 
bution is the axial marching station JX and the marching distance X(JX) . After 
the temperature distribution are two groupings which describe the extent of 
the supersonic and subsonic domains. In the first grouping the IJ refers to 
the z-grid location. The value of ISS is the y-grid location of the first 
supersonic grid point relative to surface 1 while ISSD is the y-grid location 
of the first supersonic grid point relative to surface 2. In the second 
grouping IJ refers to the y-grid location. ISS is the z-grid location of the 
first supersonic grid point relative to the 3 surface while ISSD is the z-grid 
location of the first supersonic grid point relative to the 4 surface. A 
value of ISS or ISSD of 999 refer to a row or column of all subsonic points 
and a value of ISS or ISSD of 1000 means that the corresponding surface value 
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o£ the Mach number is supersonic unless that surface happens to have a wall 
function boundary condition. Following the initial conditions, the above 
procedure is repeated at each marching station until the station is 

reached. At each station where a restart is written, a message is written. 

The message tells the file name, the marching station number and the sequence 
number. The only difference between the above output and that for a laminar 
case is the lack of printout associated with the turbulence related variables. 
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APPENDIX A 

THE ROBERTS' TRANSFORMATION 


An analytical coordinate transformation, developed by Roberts (Ref. 23) , 
is a very effective means of concentrating grid points in the region of the 
computational boxandaries. If, for instance, the region near a computational 
boundary corresponds to a boundary layer, the concentration of grid points 
in the physical plane will help in resolving the large gradients. Suppose 
that N grid points are used in the range - X - X 2 same applies to the 

z-direction) , and that the steep gradients are anticipated in a region 3 (Y 2 “ 
near y^. Then Roberts' transformation y,j,(y) is given by 




N + (N 



b+y-c 

b-y+c 



b+a ^ 
b-a / 


(A-l) 


2 2 

where a. = - y^^, b = a /(1-3) and c = y^. The use of equally-spaced points 

in the transformed coordinate, y^, ensures a concentration in the region 

S. y<Yj^'*’3(Y2 “ • Derivatives with respect to the physical coordinate, y, 

are obtained from the chain rule formulae; 


and 


d d 

ay “ dy ay^ 


(A-2) 


ay^ " \ ay j "ay^ ay^ dy^. 


The relationships for dy^/dy and d y^/dy 


are easily calculated from Eq. (A-l). 
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APPENDIX B 

SPECIAL BOUNDARY CONDITIONS 

Most of the boundary conditions utilized in the PEPSIS computer code are 
fairly straightforward and have been included in the users ^ manual portion of 
this report without comment. Three types of boundary conditions not in 
common usage are also available to the PEPSIS user. They are (1) Mach line 
extrapolation, (2) wall functions and (3) the momentum equation as a boundary 
condition. 

Mach line extrapolation is predicated on the concept that in a simple wave 
region, the flow properties remain constant along Mach lines. For a scalar 
quantity, (|), this results in the equation 

*>* 

while for the velocity, V, this results in the equation 

n^-7v-0 

where n is the unit vector in the direction of the Mach line. In two 
m 

dimensions, the Mach angle relative to the coordinate system, ip, is given by 
the sum of the flow angle relative to the coordinate system and the Mach angle 
relative to the flow; viz. 


t-tan-'(-^) + sin-'(^) 


Thus in general orthogonal coordinates Eq. (B-1) becomes 


d<f> 

dx“ 


h , d(f> 

+ + =0 


while Eq. (B-2) , a vector equation, yields two scalar equations 


(B-3) 


(B-4) 
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du 

It 


(B-5) 


+ 



dh, h, / au 

-r^v + -J- tanU/f-r — 

ay h ay 



- 0 


and 


dv 1 ah, h, / av I ah, x 


(B-6) 


The Mach line extrapolation boundary conditions Eq. (B-4) through (B-6) are 
typically applied on the boundary upstream of the cowl lip. Eq. (B-4) is 
associated with the continuity equation (with <j> = P) while Eq. (B— 5) and 
Eq. (B-6) are associated with the x- and y-momentum equations, respectively. 

In the PEPSIS computer code, wall functions may be used for the wall 
boundary conditions in lieu of the no-slip condition. Wall functions (as 
utilized in PEPSIS) are valid only for turbulent flow and are typically 
used when it is not possible to have enough grid points in the wall boundary 
layers to adequately resolve the large gradients normal to the wall. The 
wall function concept assumes that the velocity parallel to a wall, u, obeys 
the logarithmic law of the wall profile 


U 



Iny*^ \ 

— 


(B-7) 


where 

y 

^ lam 


(B-8) 


u^ is the wall shear velocity and “y is the normal distance from a wall. Taking 
the derivative of Eq. (B-7) with respect to, say y (for a 1 surface), yields 

du (b-9) 

ay “ 
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For a 1 surface, is given by the relationship 



(B-10) 


and by assuming that the ratio of u/w does not vary (as a function of y) 
in the region of a wall, the relationships for 3u/3y and 3w/3y can be derived. 
They are 


du 




(B-ll) 


and 


w 


(B-12) 


The boundary conditions represented by Eqs . (B-ll) and (B-12) are applied at 
the first grid point away from a wall rather than at a wall grid point since 
the logarithm law of the wall is not valid at the wall. Solution of Eqs. 
(B-ll) and (B-12) in conjunction with the governing partial differential 
equations results in nonzero wall slip values for u and w. These slip 
values represent the values of u and w required to satisfy the derivative 
conditions, Eqs. (B-ll) and (B-12). Equations analogous to Eqs. (B-ll) and 
(B-12) can easily be derived for surfaces 2, 3 and 4. 

The last category of boundary condition considered in this appendix is 
the normal momentum boundary condition. This boundary condition is typically 
applied as one of the wall boundary conditions (the other wall boundary 
conditions, for instance, might be u = v = w = 0) . If a 1 or 2 surface 
is considered, Eq. (4) can be utilized, except that in this case several of 
the terms can be dropped. In the zero bleed case (v = 0) , Eq. (4) becomes 
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(B-13) 


75- ( h,h,/.v*) - h,-5p- /,u* - h, -5^ pw* 

dp 4 d / >^ 1^3 dv \ 

**■ ^1^3 dy " ” 3 ” dy \ hg ^ dy / 


For the case of u, v, w and their derivatives small, Eq. (B-13) becomes the 
boundary layer approximation 

dp ^ (B-14) 

-7 — ■ 0 

dy 

An equation analogous to Eq. (B-13) can easily be derived for a 3 or 4 surface 
except that in this case Eq. (5) is used as the starting point. 


*? 
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Figure 6. - Variation of pressure coefficient at selected spanwise locations. 
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Figure 10. ~ Cross-section of NASA Ames inlet. 
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Figure 11. - The computational mesh. 



Figure 12. - Inviscid shock pattern. 
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Figure 13. - Mach number profiles in center plane of inlet. 


















TABLE I. ~ PEPSIS SAMPLE INPUT 


NASA - LEWIS TEST CASE 

021.0 


SLIST2 


IWR = Of 


USCALE = 

1 .Of 

TZERO = 

A 1 6 • 0 f 

UZERO = 

3000. Of 

YZERO = 

1 .16666667 

VISCOS = 

0.01 f 

SEND 


SLIST3 


ICAR = 0 

f 

lAXl = 0 

f 

SEND 


SLISTA 


XJXTW = 

3.72Ab3f 

AMDOTD = 

16000*0.0 

YS ( 1 f 1 ) 

* 

o 

• 

o 

II 

YS(2f 1) 

= 1 .Of 

YSn f 2) 

o 

. 

o 

II 

YS (2f 2) 

in 

. 

o 

II 

NE(1) = 

20. 

NE(2) = 

h • 

IGRIO ( ] > 

= ?f 

IGRIU (2) 

= 0, 

EPS(l) = 

0. 1 Of 

XENTR = 

0.91B21 » 

XI - 0.0 

f 

X2 - 0 . 025 f 


70 



AP = 1.0» 

NS = 3* 

IRSTOT = 3» 

IRSTIN = 0, 

NFILE=0* 

NSAVEO=Of 

5END 

SLIST5 

ILAM = 0* 

SEND 

SLIST6 

YSL0T(2) = 0.05» 

SEND 

SLI5T7 


NRCON ( 1 

• 1 

» 1 ) =4, 

N3C0N ( 1 

t2 

II 

N3C0N ( 1 

»3 

» 1) =1 » 

NBCON ( 1 

*4 

^ 1 ) =1 * 

NBCON (2 

f 1 

» i ) =0f 

NBCON (2 

»2 

f 1)=1» 

NBCON (2 


» 1 ) =1 » 

NBCON (2 

.4 

, 1 ) =1 » 

NBCON ( 3 

♦ 1 

f 1 ) =1 » 

NBCON (3 

♦ 2 

» 1 ) =1 » 

NBCON ( 3 

»3 

»1)=1* 

NRCON (3 

*4 

» 1 ) = 1 » 

NBCON (4 

*1 

II 

NBCON (4 

f 2 

» 1 ) = 1 » 

NBCON (4 

»3 

, 1) =0, 

NBCON (4 

*4 

» 1 > =1 * 


SEND 

SLISTQ 

SEND 



“vl 

ro 


TABLE II. - PEPSIS SAMPLE OUTPUT 





DiFre^ENCE: 

DPE^ATQf^S 

OIRECTION -2 



I 

Y(I ) 

AG (I, 1 ) 

AG{ 1,2) 

A G C 1 , 3 ) 

aSH,4) 

AG tlt5 ) 

AG(lt5 > 

1 

0. 0000 

-2.9429+01 

1 . 1772+02 

- 8 . 8288+0 1 

3.9352+03 

- 8 . 8 122*03 

4.8770*03 

2 

1.9489-02 

2.2538+01 

L.3000 

-2.2538+01 

1 .766 q+03 

-4.0638*03 

2.2973*03 

3 

4.4783-02 

1 .7471+01 

0.30Q0 

-1,7471+01 

1 .3699+03 

-2.4420*03 

1 .3721*03 

4 

7.7159-02 

1 . 57 + 01 

0,3000 

-1.3757+01 

6 .6997+02 

-1.5 140*03 

8.4399*02 

5 

1.1787-01 

1 .1349+01 

3. 3 0 00 

-1.1049+01 

4 .3756*32 

-9,7656* 02 

5.3900*02 

6 

1 . 6 794-01 

9 ,0949+00 

0.0000 

-9 , 0949+ aQ 

,3.ai3Q.*02_ 

- 6 . 6 174*02 

3.6074*02 

7 

2. 2787-Gl 

7.7139+DO 

Q.OOGD 

-7.71 39+00 

2 .204 1 + 02 

-4. 7604*02 

2 ,5563^02 

8 

2.9729-01 

5 , 7770+00 

3.0000 

-5.7770*0 j 

1 .7359 + 02 

-3.6743*02 

1.9384*02 

9 

3. 7474-01 

5.1970+00 

«J • Lt G w Ij 

- 6 . 1970 + no 

1 .4838+32 

-3. 0722*02 

1.5884*02 

10 

4.5761-01 

5.9199+03 

c.oooo 

-5,9199+00 

1 ,3857 + 02___ 

-2,8 036*02 ■ 

1.4180*02 

li 

5.4259-01 

5.9199+00 

0.0003 

-5.9199+ OD 

1 .4 18 0*02 

-2,8036*02 

1.3357*02 

12 

6. 2526-01 

5 .1970+00 

30.00 . . 

-5*1970+0Q 

1 .5864 + 0? . 

-3.0722*02 

1.4838*02 

13 

7.0271-Ql 

5.7773+00 

C .00 00 

-6.7770+0 J 

1 .9 364+02 

-3.6743*02 

1.7359*02 

14 

7.7213-01 

7.7139+00 

0. 3 0 00 

-7.7139+03 

2.5563+02 

-4. 7604*02 

2 . 2041*02 

15 

3.32 - 0 i 

9 . 3949 + 00 

0.0000 

-9.3949+00 

3.6074+02 

- 6.6 174*02 

3.0100*02 

16 

3.821 3-01 

1 . 1 0 4 9 + Q 1 

0.0000 

-iiLlJ.49*_31 

..._.S,i9_00*0 2_.__ 


4.3756*02 

17 

9. 2234-01 

1.3757+01 

D . J 3 00 

-1.3757+01 

8 .4399+02 

-1.5 140*03 

6.6997*02 

18 

9.5522-01 

1.7471+01 

0*00 00 


1 * 3.72 1_+D_3 

-2.442D+03 

1.0699*03 

19 

9.8051-01 

2 .253S + D1 

ij.oaoQ 

-2.25 38 + 01 

2 .2973 + 03 

-4.0638*03 

1.7665*03 

?Q 

1.0003+00 

8 ,S2d8 +01 

- 1 . 1772 + 02 

2.,9 4z?L»-J)X_ 

_ 4.8770*33 

-3.8 122*03 

3.9352*03 


DIFFERENCE OPERATORS DIRE:TI0N-3 


Y( I) 

Q* QOGD 
1 . 30 . 00-31 
2 * 300 Q-Qi 
3 . OOQD-Ql 
M. 3030-01 
5 * 3030-01 


AG(I,1 ) 

'5 • OQQ J + 00 
5 . 3JKiO:*.?D.._.. 
5 . 3003+00 
S.OOOGfGD 
5 . 0OlJ3 + QO 
1 .5300 + 01 


A 0 ( I r 2 > 

2 . 00 00+ 31 

0.0030 
3.3000 
0.3003 
■ 2,3300 + 01 


A 3 C 1 , 3 ) 

1 . 5300+01 

5 

5 .OODO+Ou 
5..oaoo*oo 
5 . 3300+00 
5 .33ijo*04J 


AS (I, 4) 


AG (1,5 ) 


l.OOQD +02 - 2 . 0000+02 

1 . 300 n+G 2 - 2 . 0000+02 

1 .0 oa 0+ J12 -2._a DQIL+ Q 2- 

1 . 3300 + 02 - 2 , 0000*02 

1 , 0Qail +_ Q2 OaPD * Q2 . 


aS*I»5> 

1.3000+02 


1.0000+02 


1 . 0000+02 


'vl 

00 










METRIC 

DAT4 





1 

X 

Y 

FG INHX > 

F0(NH2I 

FG(Nh3> 

F3t nH12 ) 

FGtNH21 » 

FS (MHM) 

F6|NH32) 

BETA 

20 

9,182-01 

1.000*00 

9.831-01 

9.931-01 

1. 000*00 

-5*326-04 

-3.628.rQ2_ 

_0».QOO 

Q^.aaa 3.664-06 

19 

9.182-01 

9.805-01 

9.831-01 

9.331-01 

1 .300*00 

-l.b2r-03 

-3.626-02 

□ . 000 

0.000 

4.793-02 

IB 

9,162-01 

9. 552-01 

9.832-01 

9. 932-01 

1.000*00 

-3. 726-03 

-3.619-02 

0.030 

Q.QOO 

1 . 101-01 

17 

9. 1 82-01 

9.228-01 

9.833-01 

9.333-01 

1 .000*30 

-6.441-03 

-3.604-02 

0.000 

0.000 

1.B9S-Q1 

16 

9.182-01 

P.621-01 

9. £37-01 

9.337-01 

1 .000*30 

-9'. 3 89-03 

-3.5h-i)2: 

__0. OOQ 

_fltQaiL 


15 

9.182-01 

8,321-01 

9,843-01 

9 .9 4 3-01 

1 .000*30 

-1.420-02 

-3,512-02 

0.030 

o*oOO 

4,399-01 

m 

9.182-01 

7. 721-01 

9.853-01 

9.353-01 

1 .OOJ 4 OO 

- 1 . 954-02 


_ _Q. DOd 

O.QOQ 

13 

9. 182-DI 

7.027-01 

9.869-01 

9.359-31 

1 .300*00 

-2.604-02 

-3.243-02 

0.030 

0.000 

T,l3i-01 

12 

9.182-01 

6.253-01 

9.692-01 


i .000*30 

t3, 367-Q2 

*9S5-’02 

□ . 00 Q 

n.nno 

8.844-01 

11 

9.162-01 

5.424-01 

9.924-01 

9.924-0 1 

1 .000*30 

-4.314-02 

-2.603-02 

Q.ood 

0.000 

l.Q5S«00 

10 

9.182-01 

4.576-01 

9.965-01 

9.965-01 

1 .000*30 

-5,397-02 

-2 *061-02. 

_. 0*_000 ... 

D *000 

l.glQ9DQ 

9 

9.182-01 

3.747-01 

1 .oci* oc 

1.301*30 

1 , 000*30 

-6.633-02 

-1 .329-02 

0*030 

0.000 

1.337*00 

8 

9. 182-qI 

2.973-01 

1 . ^07*00 

1 .oo7«oa 

1 .000*30 

-e. 001-02 

-3.980-03 

,__D*DOO 

Q»OOQ- 

1. 4259(10 

7 

9. 182-Ul 

2.279-01 

1 .U13* CO 

1.013*00 

1 .303*00 

-9.459-02 

7,144-03 

'0. 000 

0.000 

1.469*00 

b 

9.182-01 

1. 679-oi 

1 . C 1 9* OC 

1.319*30 1,000*30 

-1.394-01 

1 .9 i;u-D2 

0.000 

O.OQD 

_1.4TS*0fl_ 

5 

9, l8t-0l 

1.179-01 

i .C25*0C 

1.325*00 

1.300*00 

-1.258-01 

3,249-02 

0. GOO 

Q.OOO 

1 . 449*00 

4 

9. 182-01 

7. 716-32 

1 .G30*QG 

1.330*30 

1 .300*00 

-1.371-01 

4.5l8ro2 

0*000 . 

Q.OOD 

„_U4a4*0Q. 

3 

9.132-al 

8. 470-02 

1.036*00 

1 ,335*30 

1 .000*30 

-1.492-01 

5.710-02 

0. 000 

0.000 

1.34 7*00 

2 

9. 182-01 

1.949-02 

1 .D39*GC 

1 .339*00 

1 .000*00 

-1.550-01 

. 6. 906-02 

0.000 

0.000 

U2A6U10. 

1 

9.182-01 

C . 0 0 Q 

1.04 1*03 

1.041*30 

1 .300*00 

-1. 148-01 

7.972-02 

0. 030 

0.000 

1.086*00 



JXr 

1, VARIABLE 

"4 

X 

1 1 

.91B21+DC 





.... 



LEVEL 3 

♦ ♦ ♦ * ♦ 

U-VEJ^ _ 

***** 



iz - 

1 

2 

3 

4 

5 

6 


2 z 

.3000 

. luQO+03 

•2000+00 

. 3Q3Q+0D 

.4000+00 

.5000+00 

20 

•1D03+01 

.1000+31 

. lQCO+01 

. 10D3+0 1 

. lDDO+01 

.1000+01 

• IDOO+Ol 

19 

•93 13+DG 

.1000+01 

. It. 00 + ul 

.1000+01 

. 1000 + 0 1 

.1000+01 

. 1003+31 

18 

•9561+00 

.1030+31 

. ICCO+Ol 

.1000+31 

. lCDO+31 

.1000+01 

. 1000+01 

17 

•9243+00 

.1000+01 

. 1GD3+01 

,1303+01 

. 1000+01 

1 0.0 0+^1. _ 

.1030+01 

1& 

•3842+00 

.1003+31 

.1030+01 

.1300+01 

. 1300 + 01 

.1000+01 

. 1030 + 01 

15 

.3349+00 

.1003+31 

. 1 0C0+ 01 

.1300+01 

. 1300+01 

.1000+01 

, .. 1000^01 

14 

,7 7 5.9 + 00 

.1030+01 

. 13D0 + O1 

, 1000+01 

• looo+oi 

.1000+01 

. 1000+01 

13 

.7075+00 

.9999+00 

. 9999 + 3C 

. 9799+00 

. 'J9 99+Q0 

.9999+00 

• 9999+On 

12 

.5310+00 

.9999+00 

. 9999+ 03 

.9999+03 

.9999+00 

,9999+00 

.9999*00 

11 

.5489+OC 

.9998+30 

.9993+03 

,9?98+00 

.9998+00 

,9998+00 

.9998+00 

ID 

.4645+03 

,9998+00 

.9998+03 

.9998+00 

. 9998 +03 

.9998+00 

• 9998 + 00 

9 

.3818+00 

.9997+3L 

.9997+03 

.9997+00 

,9997 +DC 

.9997+00 

.9997+00 

8 

. 3040+00 

.9997+30 

.9997+03 

.9997+00 

. 9997+03 

.9997+00 

.9997+00 

7 

.2339+00 

,9997+3 J 

, 0997+ QC 

. 9997+30 

. 99 97+3 0 

.9997+00 

.9997+00 

6 

• 1 7 3u + L»Q 

.9997+00 

. 9997+CD 

,9997+30 

, 999 7 + 00 

,9997+00 

.9997+00 

5 

. 1218+OD 

.9997+30 

. 9 997+ CO 

,9997+OD 

.9997+30 

.9997+00 

.9997+00 

4 

.7993-01 

.9997+00 

.9997+30 

.9997+33 

,9997+33 

.9997+00 

.9997+00 

3 

.465C/-C1 

,9997+00 

,9997+30 

.9997+30 

. 9997+30 

.9997+00 _ 

. 9997+DO 

2 

.2027-01 

.9993+00 

.9998+00 

.9998+30 

.9993+00 

,9998+00 

.9998*00 

1 

.3000 

.9998+00 

.9993+00 

, 999B+0^„ 

, 99_9_8 +.0 0_ _ 

.9998+00 

.9998+00 


•vl 

(J1 



■vl 


JX r 

1 , variable: A 1 XI i» 

= .91321+00 







LEVEL 

3 tc 4 c 4 c ♦ « 

V-VEL 




IZz 

1 

2 

3 

4 

5 

6 


z = 

.0030 

. luD3+00 

.233 0 + q 0 

• .?aoo_io i3_. 

.4300+00 

.5000+00 

2 Q 

.iDOD+ai 

-.6394-07 

-.6394-07 

-. 6594-37 

-. 63 94 -3 7 

-.5394-07 

-.6394-07 

19 

.9BiO+ua 

-.3366-33 

-.8266-03 

-.8366-33 

-.83_66-03 _ 

-^3366-07 

— .8366— 03 

IB 

.95b1-»00 

- .192 1-02 

-.1921-02 

-. 192 l-o2 

-.1921-02 

-.1921-02 

1921-02 

17 

.9243+DC 

- .3307-02 

-.33Q7-Q2_ 

-.3307-32 

-.3307-32 

-.3307-02 

-.33o7-o2 

16 

.8842+00 

- .504 1-32 

-. 5041-02 

-.5341-32 

- • 53 4 i -Q 9 

-.5'34 1-02 

-.5041-02 

15 

.8349+00 

-.7154-02 

-.7154-02 

-.7154-02 

-.7154-02 

-.7154-02 

-.7154-02 

14 

.7759+00 

-.9642-32 

-. 9642-02 

-.9642-32 

-.9642-0? 

-.9642-02 

-.9642-02 

13 

.7075+00 

-.1245-01 

-.1245-01 

-.1245-31 

-.1245-01 - 

■ - iJ 2_4-5’0 1_ 

-.1245-01 

12 

.5510+00 

- . 1 54 3-01 

-. 1543-01 

-.1543-31 

-. 1543-31 

-.1543-01 

1543-01 

11 

.5489+00 

-.l84i-0i 

-.1841-01 

-.1341-31 

184 1-31 

-.1841-nl 

-. 1841-01 

iQ 

.4645+00 

-.2112-01 

-.2112-01 

- . 21 12-01 

21 12-01 

-.2112-01 

2112-01 

9 

.3818+00 

-.2333-01 

— *2 333 — ol 

-.2333-31 

-. 2333-01 

-.2333-01 

-.2333-01 

8 

.3j4q+o3 

“. 24 S 6~3 1 

-.2486-01 

-.2486-31 

-.2486-01 

-.2486-01 

-.2486-01 

7 

.2339+00 

-.2564-01 

-.2564-01 

-.2564-01 

-.2564 -D 1 


-.2564-01 

6 

.1730+00 

-.2574-01 

-. 2574-01 

-.2574-01 

-.2574-01 

-.2574_01 

-. 2 ^ 74-01 

5 

. 1 218+ JO 

- .2529-01 

- , 2529-01 

-.2529-31 

2529-11 

-.2529-01 

-. 2529-01 

4 

.7993-01 

-.2450-01 

-,?453-oi 

-.2450-31 

-.2450-01 

-.2450-01 

-.2450-01 

3 

.465Q-Q1 

-.2353-01 

— • 2353 — ul 

-.2350-31 

-.23 50-31 

-.7350-01 

-.2350-01 

2 

.2027-01 

-.2209-01 

-.2209-01 

-.22q9-o1 

-. 22 09-01 

-.2209-01 

2209-01 

1 

.0300 

-.1895-01 

-.1895-01 

-.1395-31 

-.1895-01 

-.1395-01 

-.1895-01 



JXZ 

It variable 

AT XC 

1 )- .91821 + 03 

- - - --- 









level 3 


w-veT " 




U- 

1 

2 

5 

4 

5 

6 


Z r 

.0030 

, 1000+ DO 

.2330+30 

.30 DO +30 

.4300+00 

.5000+00 

20 

• 1 000+Qi 

.3000 

.Ol^QD 

,03 00 

. 00 30 

__ «0.3Q0 

.0000 

19 

•981u+00 

.3000 

.0030 

.0330 

.0000 

.3300 

.0000 

18 

. 9561*00 

• 0030 

.0003 

.0330 

.GOOD 

*3.inn 

.0000 

17 

.9243+00 

.3000 

. 0003 

• 0300 

. 0000 

.3300 

.0000 

16 

•3842+00 

.0000 

.0000 

.0300 

.0000 

• 3300 

• 0000 

15 

.8349+00 

.0030 

• 0 G 0 3 

.0330 

. 0000 

.3300 

• 0000 

m 

.7759+00 

•3030 

* Q u 0 0 

.0300 

.0300 


.3000 

13 

.7075+00 

.3000 

.0033 

. 3330 

. 0300 

.3000 

.0000 

12 

.6313+00 

.3030 

.0G03 

.0300 

. 030D 

•3300 

.0000 

11 

.5489+00 

.3030 

. OQDO 

,0333 

, GOOD 

.3300 

.0000 

10 

.4645+00 

• 0000 

.0300 

.0300 

.GOOD 

• 30nn 

.0000 

9 

.3818+00 

.0030 

.3003 

.3300 

. GOOD 

.3300 

• GOOD 

6 

.3340+00 

.3030 

.0003 

.0330 

. GDOO 

. 3300 

. Onnn 

7 

.2339+00 

.3000 

.0 003 

.0330 

.GOOD 

.3300 

.0000 

6 

.1730+00 

.3000 

. n c 0 3 

.0300 

.00 00 

.3000 

■ .0000 

5 

.1218+00 

.3030 

.0033 

. 0330 

.0300 

.3300 

.0000 

4 

.7993-01 

.3030 

.0003 

.0330 

•.QQ.Cm._.. 

.3000 

.0000 

3 

. 465u~0 1 

.3000 

.0300 

.0300 

.00 00 

.3300 

• OOQg 

2 

.2027-01 

.3000 

. G G 3 

.0300 

.COOQ 

.3300 

• 0000 

1 

.3000 

.3030 

.0003 

, 3300 

....DQIliL 

.3300 

.0000 


VJ 

VJ 



■vj 

00 


JX r 

If VARIABLL 

AT X ( 1 ) = 

.^1821+00 







le:ve l 

3 

RHO 





1 


3 

4 

5 

6 


Z - 

.0000 

. 1300+ GO 

. 2300+03 

. 3000+00 

.‘4000+00 

. 5000+00 

20 

. lODU+Dl 

. 1 D 0 Q + Q 1 

.1003+01 

.1330+31 

. 1300+01 

.1300+01 

.1000+01 

19 

.981c, + 0D 

.1000+01 

.1^:00 + 01 

. 1300 + 31 

. 1000+01 

.1000+01 

. 1000 + 01 

18 

•9S61 +Q0 

.1000+01 

• 1 GCO'^ Ol 

.1303+31 

. 1000+01 

.1000+01 

. loo^+m 

17 

,9243+00 

.1000+01 

. IGOO+Dl 

.1300+01 

. lOQO + jl 

.1000+01 

. 1000+01 

16 

.3842+00 

.1030+01 

.1300+01 

. no 0+01 

. 1000+01 

.1300+01 

. 1000+01 

IS 

.8349+00 

•1 OOO+Ol 

. 1 jOD+ 01 

. 1300+0 1 

. 1000+01 

.lOOc+Dl 

. 1000+01 

1** 

*7759+00 

.1030+01 

. K QO + Ol 

. iOOO + Ol 


.1300+01 

.1000+01 

13 

.7075+00 

.1030 + 0 1 

. KG3+Q1 

. 10 D 0 + 3 1 

. 1000+31 

.1000+01 

. 1000+01 

12 

.6313+DO 

.1000+01 

. 1 JQ3+ Ql 

. 1300 + 31 

. . 1DDQ.+31_ _ 

.1000+01 

.J-0-DQ*ai 

11 

,5469+00 

. 1030+01 

. IDOO + ul 

. 130D+3U 

. 1000 + 3I 

.1300+01 

. 1000+01 

10 

.4645+00 

.1030+31 

. 1000 + 01 

. 1300+dl 

. lODO+Ql 

.1300+01 

. 1000+01 

9 

.3818+00 

.1 DGO+01 

. IGQQ+ai 

. 1300 + 3 1 

. 1000+01 

.1 300*(31 

. 1000+01 

B 

.3040+00 

.1030+01 

. 1000 + 01 

. 130D+nl 

. lODQ + 01 

.1300+01 

• lOOO+Ol 

7 

•2339+uO 

.1000+31 

. 1000+ Q1 

.1300+04 

. lQOQ + 01 

.1300+01 

• lOOQ^Dl 

6 

. 1 7 3 rj + Q 0 

.1000+01 

. 1000+01 

. 1300+01 

_ 

.1300+01 

.1000+01 

5 

.1218+00 

.1003+01 

. 1 ODD'* d 

.1300+01 

.1000 +01 

.1300+01 

• 100 Q+Ql 

4 

.7993-01 

.1000+31 

. 1000 + 01 

.1300+31 

. lDOD+01 

.1 ODD+Ol 

. 1000+01 

3 

.4650-01 

.1030+01 

. 1 000 + 31 

.1300+01 

. 1000+01 

.1300+01 

. 1000+01 

2 

.2027-01 

.1 QjQ+Ol 

. 1 000+01 

.1303+01 

?.iO_LQ_+_3 l-„„ 

■ la-OB^ai 

.1000+01 

1 

. 0 0 0 ij 

. 1 030+01 

. IL’ 00+01 

. 1330+r)l 

. lOOQ+Ol 

.1300+01 

• 1000 + 01 



JX r 1 , V A R I A b L E ' AT X ( i ) = .91 8 2 1 + 0 6 






ir^EL 

3 

ENERSY 




IZ- 

1 



3 

4 

5 

6 


Z r 

.3030 


. 1033 + 03 

.2330+30 

. 3000+30 

.4300+00 

• 5000400 

20 

.1300401 

. 3 0 3 Q 


.0003 

.Q3D0 

.3000 

.3300 

• 0000 

19 

.?aiu+QC 

.3030 


n . r» T 

t vJ C. ij iJ 

.0303 

. 3000 

.3300 

.0000 

18 

•9561+uQ 

.3030 


.0003 

.0300 

. CDOO 

.0300 

.0003 

17 

.9243+00 

.3D3l3 


.0003 

. 03 30 

. 0000 

.3300 

.0000 

16 

. B842 + C0 

.0030 


.3300 

. 3300 

.0030 

.3300 

• 0000 

15 

.3349400 

.3000 


.3000 

.0300 

.0300 

.3330 

.0000 

m 

.7759400 

.3030 


. Q 0 0 3 

.0300 

.CODQ 

.3300 

.0300 

13 

.7075+00 

.3030 


.3000 

. u 3 0 0 

. ODQD 

.3300 

• OQQQ 

12 

.6310+00 

.3030 


.3003 

. 0300 

.3000 

.3300 

*0000 

11 

.5489400 

.3030 


.3 300 

.0300 

■ .00.3u 

*3300 

•OOQD 

10 

.4645+00 

. j 0 3 0 


.0003 

. 0303 

. uG 0 U 

.3300 

.0000 

9 

.3B1&400 

.3030 


.3033 

.0300 

-0300 

.3300 

.0030 

8 

.3040+03 

.3030 


. 3CU0 

.0303 

.0300 

• 3300 

• OOOO 

7 

.2339400 

.3030 


.3003 

.0333 

. 0003 

.3300 

. 0000 

6 

.1730400 

.3030 


. 3 0 3 

.0333 

. aao-Q- 

. . innn 

-..■DO on 

5 

. 1 218+uD 

•3030 


. 33 GO 

.0303 

. 03 00 

• 33qo 

.0000 

4 

.7993-01 

. 3030 


.3000 

. 3303 

. oano 

-3300 

- . 00 an . 

3 

.4650-01 

.3030 


.0033 

.0333 

. 00 00 

.3300 

.OOgo 

2 

. 2 0 2 7 — u 1 

.3300 


. Q 0 3 0. 





.3300 

■-.Qpnn 

1 

. 3000 

.3000 


. 3033 

.0300 

. 0330 

.3300 

.0000 


•Vj 

VO 



00 

o 


JX- 

1, VftRIABLt: 

A T XI 1 ) = 

.91821+00 








LEVEL 5 






12- 

1 


5 

4 

5 

6 


z - 

.3000 

. 1030 + 00 

*^’30 0 + 3 0._ 

. » 3jjD+Q.O 

.4300+00 

. *505.0^0X1 

2G 

• nnj + Gi 

. 1 D 3 0 + 1 1 

. 1000+ D1 

. 1303+01 

. 103Q+01 

.1300+01 

•l000*0l 

19 

. 9 B 1 j + J 0 

.1003+01 

. ILOO+Ol 

. 1330+31 

• lOOp+Qi 

. 1 nnn+Qi 

-+.lD.Qil+Q.l 

le 

. 9561 + JO 

.1300+01 

. 1 J03+D1 

• i330+3l 

. lQOO + 0 1 

.1000+01 

. 1000+01 

17 

. 9243+ QQ 

,1003+31 

. laoo+Gi 

.1300 + 3 1 . 

. lQOD+31 

- .logn+Dl- 

.IDDQ+Ql 

16 

. 6842 + 0 'J 

.1003+01 

. i .j'no + oi 

. 1D3Q+31 

. 1000+31 

.1300+01 

. 1000+01 

15 

.8349+00 

.1333+01 

. 13C3+ui 

. 1300+01 

. 1030+01 

.1300+01 

. 1000+01 

14 

.7759+00 

.1033+31 

.1300+01 

.1330+31 

. lODG + 01 

.1300+01 

. 1003 + 01 

13 

.7D75+LQ 

.1033+31 

• 1 C 3 0 ♦ 0 1 

. 1333+31 

. 1300+01 

.1300+01 

. lOno+oi 

12 

. 8 3 1 0 + 0 0 

. 1030+01 

. 1003 + 31 

. 1300+31 

. 1000 + 31 

• 1300+01 

. 1000 + 01 

11 

.5489+00 

.1000+01 

.1000+01 

.13Q0*31 

. 1300+01 

.1300+01 

• 1000+01 

10 

.4645+00 

• 1 DJD + 3i 

. 1300+ 01 

. 1303+31 

.1000+31 

. 1 OOo+Ol 

.1000+01 

9 

. 3 BI 8 +OC 

.1003+31 

• 1 GOO + 01 

. 1300+31 

.1300+31 

.1300+01 

• lDDQ+01 

8 

.3040+00 

.ICDD+Oi 

. 1003+01 

.1300+01 

. 1300 + 01 

.1300+01 

• IDoo + Ol 

7 

.2339+30 

.1 330+01 

. 130D+D1 

.1330+31 

. lopD+ai 

.1300+01 

. 1000 + 01 

6 

.1733+00 

.1033+01 

. 1003+ Ql 

.1333+31 

. 1000+01 

.1300+01 

•1000+01 

5 

.1218+30 

.1000+31 

. 1000+ 01 

. 1330+31 

. 1000+01 

.1300+01 

. 1000+01 

4 

. 7995-31 

• 1 030+01 

.1003+01 

.1330+31 

. 1300+31 

.1300+01 

.1000+01 

3 

.4650-01 

.1003+01 

. 1003+01 

• llOQ+Ql 

. IQOO+Ol 

.1000+01 

•1000+01 

2 

.2327-01 

.1003+01 

.1000+01 

. 1300+31 

. 1000+31 

•laoo+oi 

.1000+01 

1 

.0303 

.1000+01 

. icno + Di 

•i^oatoi 

.iaoo^+jii„ 

.1300+01 

.1000+01 



JXr 1, VARIABLE AT X( 1>= .91821+00 

LEVEL 3 



iz- 

1 

O 

t- 


Z - 

.0000 

. 1 G 00 + on 

20 

.1000+01 

.3122+00 

.3122+00 

19 

.98 lo+on 

.3122+30 

.3122+00 

18 

.9561+DQ 

.3111+30 

.3111+00 

17 

.9243+ao 

.3093+00 

.3093+00 

16 

.5842+00 

.3063+30 

. 3C65+CD 

15 

.8349+00 

.3q1 3+OD 

.3013+00 

m 

.7 7 5^ + 00 

.2936+00 

.2936+00 

13 

.7075+00 

.2820+30 

.2320+00 

12 

.6310+00 

.2648+00 

.2643+ OC 

11 

.5489+00 

.2436+03 

,2406+00 

10 

.4645+00 

.2074+00 

. 2074 + 00 

9 

.3818+00 

.1640+00 

.1643+00 

6 

.3040+00 

.1101 +30 

. 1 101 + 00 

7 

.2339+00 

•471B-Q1.. 

.4718-01 

6 

.l7 30+00 

.2223-01 

. 2223-01 

5 

.1218+00 

,9468-01 

.9468-01 . 

4 

.7993>ol 

.1693+00 

. 1693+ 00 

3 

.4650-01 

.2879+00 

.2879+00 

2 

.2027-di 

,6 445 + Qu 

. 6445 + 00 

1 

.3000 

-.1178+01 

_ n 7_8 +_Dl__ 


00 

H 



00 

ro 


JXr 

i, variable 

AT X t i}- 

,91 8 21+ DO 

• - - 

— 






LEVEL 3 


TKe 





1 

u. 

3 

4 

5 

6 


z = 

.GOOD 

. 1000+ CO 

.2000+00 

• 3QnQ-»“0D 

. 4300+00 

,5000400 

20 

• 1 3DQ + 1 

.1501-G3 

. 1501-03 

.1501-03 

. 1501-33 

.1501-03 

.1501-03 

19 

.981 J+00 

.1531-03 

. 1 501-03 

. 1501-03 

, 1501-03 

.1501-03 

. 1501-03 

18 

.9561+UD 

.1531-03 

. 15Q1-C3 

. 1501-03 

. 1501-03 

.1501-03 

. 1501-03 

17 

.9243+QC 

, 1531-33 

.1501-03 

. 1501-03 

.15 01-03 

.15q1-o3 

. 1501-03 

1^> 

.56i42*0D 

.1501-03 

. 1 501-03 

.1501-03 

, 1501-03 

,1501-03 

. 1501-03 

15 

.33M9+CC 

.1501-03 

.1-01-03 

,1501-03 

. 1501-03 

.1501-03 

• 15[jl-o3 

m 

.7759+00 

.1531-03 

.1501-03 

.1501-03 

. 1501-03 

.1501-03 

. 1501-03 

13 

.7D75+DC 

.1501-03 

. 1501-03 

,1501-03 

. 1501-33 

.1501-03 

. 1501-03 

12 

.&31Q + iJ0 

. 1 5 3 1 -0 3 

. 15D1 -QZ 

.1501-03 

.15 01-03 

.15nl-n3 

. 1501-03 

11 

.5489+cn 

.1501-03 

. 1501 -03 

. 15 ol -03 

. 1501-03 

.1501-03 

. 1501-03 

10 

.4645+00 

.1501-03 

. 1501-03 

. 1501-03 

. 1501-03 

.1501-03 


9 

,38 iB+uO 

.1500-03 

. 1500-03 

. 1500-03 

• 15DD-Q3 

,1500-03 

- 1500-03 

8 

.3040+00 

.1530-03 

. 1500-n3 

.1500-03 

. 1500-03 

.1500-03 

. 1500-03 

7 

.2339+00 

.1500-03 

. 1 500-03 

. 1500-03 

. 15Q0-Q3 

.1500-03 

. 1500-03 

6 

. 1730+00 

.1500-03 

,1500-03 

.1500-a3 

. 1500-03 

.1500-03 

. 1500-03 

5 

.1218+00 

.1500-03 

.1500-03 

. 1500-03 

. 15 00-0 3 

,1500-03 

• 1500*03 

4 

,7993-01 

.1500-05 

, 1500-03 

. 1500-03 

_ viSno-Di^ 

.1500-03 

, 1500-03 

3 

, 46 5u “C 1 

.1483-03 

.1483-03 

.1453-03 

. 1483-03 

.1483-03 

. 1483-03 

2 

.2027-01 

,5 365-04 

. 5365-04 

.5365-04 

. 5365-04 

.5365-04 

.5365-04 

1 

.3000 

.5365-04 

, 5 365-04 

.5365-04 

. 5365-04 

.5365-04 

.5365-04 


JX- 

If variable 

AT X ( 1 J- 

.913 21+00 







LEVEL 3 


DIS 




IZ- 

1 

n 

3 

4 

5 

6 


Z r 

.3000 

. 1000 + uO 

.2003+30 

• 3300 + 00 

.4300+30 

.5D0D+QQ 


.IDDD+Gl 

•o717-q4 

.6717-04 

.6717-34 

.6717-04 

.6717-04 

.6717-04 

19 

.98 1 J+DQ 

.6717-04 

.6717-04 

.6717-04 

• §717-0 4 

. ■ .6 Ti-XrP-4 

.6717-04 

18 

.9561+00 

.5717-04 

.6717-04 

.6717-34 

.6717-34 

.6717-04 

.6717-04 

17 

.92M3+00 

.6717-04 

.6717-04 

.6717-04 

.6717-34 

.6717-04 

.6717-04 

16 

.8842+00 

,67i7-Q4 

. 67l7-0.4_ 

j.671ItQ4 

.6.7 1.7-0 4. 

.5XIX-.^ 

.671 7-n4 

15 

.B349+C0 

.5716-04 

.6716-04 

.6716-34 

. 67 1 6 —3 4 

.6716-04 

.6716-04 

m 

.7759+00 

.6716-04 

.6716-04 

.6715^34 

. 6716-34 

■571 6-04 

,67i6-04 

13 

,7Q75+aQ 

.6715-04 

.6715-04 

.6715-34 

.6715-04 

.6715-04 

.6715-04 

12 

.6310+00 

.6715-04 

. 6715-C4 

.6715-04 

.6715-04 

,6715-04 

.6715.-0.4 

11 

.5489+00 

.6713-04 

.6713-04 

.6713-34 

. 6713-34 

.5713-04 

, 67 1 3-04 

10 

,4645+00 

.6712-04 

.6712-04 

..6 712-3 4 

.6712-04 

.5712-04 

.6712-34 

9 

.3816+00 

.5711-04 

.6711-04 

.6711-04 

.6711-04 

.5711-04 

.6711-04 

8 

.3040+00 

.67Q9-D4 

,6709-04 

.6709-04 

.6709-04 

- 6 7x39-04 — 

. 67D9~D^ 

7 

.2339+00 

.5708-04 

.6708-04 

.67q8-o4 

.6708-04 

.5708-04 

.6708-04 

6 

.1730+00 

.6 70 8-1x4 

*^6 708-04 

. 6X08-04 

.6708-04 

.5708-04 

•67n8-n4 


.1216+00 

.5709-04 

.6709-04 

.6709-04 

.67 q9 -04 

• 5709“04 

.6709-04 

4 

.7993-01 

.6711-04 

.6711-04 

.6711-04 

.6711-04. 

.5711-04 

.6711-04 

3 

.4653-01 

. 5597-04 

.6597-04 

.6597-34 

. 65 97-04 

.6597-04 

.6597-04 

2 

.2027-01 

.1496-04 

.^96-04 

,1496-04. 

.1406-04 

.1496-04 

. 1496-04 

1 

• 0000 

.2991-34 

, 2991 -j4 

.2991-34 

.2991-04 

.2991-04 

.2991-04 


00 

(jO 



c» 

4 > 


jXr i, VARIABLE X( 1 > - .91821+00 


LEVEL CP /? »»*»♦ 



IZ- 

1 

2 

3 

4 

5 



6 


Z - 

.3030 

. lGOO+00 

.2300+30 

.3000+00 

.4300+00 



5000+00 

20 

•1003+ 01 

.3000 

. 0000 

.0300 

. 0000 

•0300 



DODO 

19 

.9813+00 

-.1853-06 

-.1863-08 

-.1353-38 

- . 1863-33 

-.1363-08 

- 


1863-08 

18 

.9561+00 

-.1863-qB 

-.1863-08 

-.1353-38 

1353-08 

-.1363-08 

- 


1863-08 

17 

.9243+00 

-.1853-08 

-. 1 863-08 

-.1353-08 

-. 1863-38 

-.1 36 3-08 



1863-08 

lb 

.3842+30 

.3030 


.0300 

. GOOD 

.3000 



0000 

15 

. 3 i 4 9 + 0 Q 

-.1853-08 

-. 1663-06 

-.1553-08 

1863-08 

-.1363-08 



1863-08 

14 

. 7759 + i^n 

.3000 

.0000 

.0300 

.0300 

.3 030" 



ODOD 

13 

.7q754oo 

I8b3-0e 

1 863-08 

-.1363-08 

— • 18 6 3 - 0 6 

-.1863-08 

- 


1863-08 

12 

.5313+30 

-.1853-08 

-. 1663-08 

-.i363“08 

- * 1 8 6 3 —3 8 

-.1 8 6 3-08 

- 


1863-08 

11 

.5489+00 

-.1853-08 

-. 1865-CS 

-.1353-38 

-.1663-38 

-.136 3-08. 

- 


1863-08 

ic 

.4 6t|t,+ 30 

-.1853-38 

-.1 663-08 

-.1363-08 

-. 1863-oB 

-.1863-08 

- 


1863-08 

9 

.3818+00 

-.1863-08 

-.1863-08 

-.1363-08 

-. 1863-08 

-.1363-08 



1863-08 

8 

.3040+00 

-.1 853-08 

-.1863-08 

-.1353-08 

1863-08 

-.1863-08 

- 


1863-08 

7 

•2339+30 

-.1853-08 

-.1 663-D8 

-.1353-08 

-.1863-08 

-.1 363-08 

- 


1863-08 

6 

. 1 7 3 D + 0 D 

-.1853-38 

-.1663-08 

-.13 5 3-08 

-.1863-08 

-.1863-08 

- 


1863-08 

5 

.1216+00 

-.1853-08 

-. 1 863-Cp 

-. 1 3 5 3-08 

-. 1863-03 

-.1363-08 

- 


1863-08 

4 

.7993-01 

-.1853-38 

-.1663-08 

-.1563-08 

-.1663-08 

-.1863-08 

- 


1863-08 

3 

.46EQ-D1 

-.1853-08 

-.1863-08 

-.1353-08 

-. 1863-06 

-.1363-08 

- 


1863-08 

2 

.2327-01 

- . 1 86 3-06 

1863-08 

-.1353-08 

-. 1863-08 

-.186 3-q8 

- 


1863-08 

1 

.0300 

-.1853-06 

-.1663-08 

-. 1363-Q8 

-• I8 6 3r08_ 

-.1 363-08 

- 


1863-08 


91 8 21* no 


JXr 1, variable AT X( 1)= 





LEVEL 7 

4i 

I^ACH 




IZ = 

1 


3 

4 

5 

6 


Z X 

. 3 C J Q 

• 1 lOD+ 3C 

.2030+30 

. 30DD+0G 

. S.4 3 D.C+D..a 

_+5.aD2.+ 00 

20 

. 1 30[)+ Qi 

.3031+01 

. 3:.ID1 + 01 

.3031+0 1 

. 30D1+D1 

.3dq1+d1 

. 3001+01 

19 

.9610+00 

.3031+31 

.3:01+01 

.3331+GLI 

. 3001+01 

.3Qai+0L_. 

j.3.Dai+ILl 

18 

.9561+00 


. 3 C 3 1 + 0 i 

.3001+01 

. 3001 +01 

•3001+31 

. 3001 + 01 

l7 

.9243+00 

.3031+31 

.:'ooj + ci 

..30.Ql + 01_ _ 

. 30i3l-+Qi._ 

.3301+01 

-..JillLUQi 

16 

.8842+00 

.3031+01 

. 3001 + Cl 

.3031+01 

.3001+01 

.3301+01 

•3ool+Ol 

15 

.3349+00 

.3031+01 

. 3l.01 +C1 

.3031+01 

. 3001+31 

.3 001+^01-. . . 

^^aaoi+Qi. 

14 

.7759+GD 

.3031+31 

.3uCl+Ql 

. 3 0 D 1 ♦ 3 1 

. 3001 +01 

.3001+01 

. 3001+01 

13 

.7075+00 

.3 03 1 + Q 1 

.3001+31 

.3001+21 

. 3301+01 

.3 20.1 + 01 

-j-3gD.lt£ll 

12 

.5310+00 

.3031+31 

.3001+01 

.3301+01 

. 3001 +0 1 

.3301+01 

. 3DD1+01 

11 

.5489+CD 

.3001+01 

. 3 dqI + Cl. . 

3001 + 01,. _ 

• 30_Oi+O.L 

.3301+01 

.3D01+nl 

1C 

.4645+Oq 

.3001+01 

.3001+01 

.3301+01 

. 3001+01 

.3001+01 

. 30Q1+D1 

9 

.3818+00 

.3331+31 

. 3lC1 + 01 

.3201+01 

. 3301+01 

.5001+01 

. 3001+01 

8 

.3Q4j+00 

• 3 DO I'^Oi 

. 3 C 0 1 +■ 0 1 

. 3031*0 1 

. 3001 +01 

.30Q1+Q1 

. 3001+01 

7 

.2335+00 

.3031+01 

.3001+01 

.3001+01 

• 33d1 +31 

• 3 00 1 + 0,1. 

. 3001*01 

6 

.1730+00 

.3031+01 

.-001+01 

.3001+31 

. 3001*01 

.3001+01 

•30q1+q1 

5 

.1216+00 

.3001+31 

. 3001 + Cl 

.3001+01 

. 30.01 +_o 3 .. 

.xaai+Di 

_>.1Q0I,+JU 

4 

.7995-01 

.3031+01 

. 3001 + 01 

.3201+01 

.3001+01 

.3001+01 

.3001+01 

3 

.4650-01 

.3031+01 

. 3C01+C1 

. 3031 + 0 1 

. 3001 +31 

..-3^a-V+ Dl 

,. 3Da i+QLl 

2 

.2327-01 

.3031+01 

.3301+01 

•3031+01 

. 3301+01 

.3001+01 

.3001+01 

1 

.ooco 

. .3031+01 

• 3 00l + C 1 

.3001+01 

. 3001 +01 

.3001+01 

. ,3pni + m 


CO 

Oi 



00 

OV 


JXr i, VARIABLE AT X( i)- .91821+00 





LEVEL 3 


TSTAT 




IZ- 

1 

2 

3 

4 

5 

b 


1 - 

.0000 

. 1C00+ GO 

.2300*00 

. 3300+30 

.4300+00 

. 5000+00 

20 

•lDDD+01 

.1000+01 

. 1 GDQ + 01 

.1300+01 

. lDDO+01 

.1300+01 

.1000+01 

19 

.981u+0Q 

.1000+01 

. lOQD+Dl 

.1300+01 

. lDDD + 01 

.1300+01 

. 1000 + 01 

18 

.9&61+&D 

.1000+01 

. 1 000+ Ql 

. 1300+31 

. lDCQ+01 

.1000+01 

. lODn+Dl 

17 

.9243+00 

•lOQD+Oi 

.1000+01 

.1300+31 

. 1000+01 

.1000+01 

. lDOD+01 

lb 

.8842+00 

.lOQD+Ql 

. 1 GOQ+Gl 

.1300+3 1 

. IDOQ + Ql 

.1300+01 

. IQQO + Ql 

IB 

.6349+00 

•1 OOC+Qi 

. 1 000+ Ql 

.1330+01 

. 1000+01 

•lOoo+Dl 

. 1000 + 01 

1*^ 

. 7759+00 

.1000+01 

. 1000+01 

.1300+01 

. lQOD + 01 

.1300+01 

. 1000 + 01 

13 

.7075+00 

.1000+01 

. 1 000+ 01 

. 1300+01 

. lODG+01 

.1300+01 

• IDOo+ol 

12 

.5310+00 

.1000+01 

.1000+01 

.1300+01 

. lDOQ+01 

.1300+01 

. 1000+01 

11 

.5469+00 

.1000+01 

.1000+01 

• i3D0+0l 

. 1000+3 1 

.1000+01 

. 1000+01 

ID 

•4b45+Co 

.1000+01 

.1003+01 

.1300+31 

• 1000+01 

.1 000+Dl 

. lDQO+01 

9 

.3818+00 

.1000+01 

.1000+01 

.1300+31 

. 1000*01 

.1300+01 

. IDOO+Dl 

B 

. 304D + D0 

.1 COO+Ql 

.1000+01 

.1330+01 

. lODO+Dl 

.lOCo+Qi 

. lDDD+01 

7 

.2339+00 

.1000+01 

.1000+01 

. 1300 + 01 

.1000+31 

.1300+01 

• 1000 + 01 

6 

.173D+DD 

.1000+01 

. 1000+ Cl 

. 1300 + 01 

. 1000+31 

.1000+01 

. lOoD+Ql 

5 

.1218+00 

.1030+01 

.1000+01 

.1300+01 

. lODO+oi 

.1300+01 

. 1000+01 

H 

.7993-01 

.1030+01 

. 1000+01 

. lODO+D 1 

. IDOG+Dl 

.1000+01 

• 1000 + 31 

3 

.4b5Q-cl 

.1000+01 

• ) ODD+ 01 

.1300*01 

. 1000+31 

.1300+31 

. 1000+01 

2 

.2027-01 

.1000+01 

. ICDO+Ql 

. 1300+01 

. 1000+31 

.1300+01 

• 1000 + 01 

1 

. Q 0 C ij 

.1 DOC+Di 

• 1 0 0 0 + G 1 

.1300+01 

. ionc+01 

.1300+01 

. IDOQ^Ol 



lAOIzl 

1 J - c 

iSSrlOQG 

I s s D r 1 n an 

IADI-1 

IJ= 3 

IS S- I ODD 

iSSOrlDOO 

IADI:^1 

d" 
1 t 

r) 

hH 

XSS-IOD j 

I SSD=icOO 

lADIzl 

Ijz 5 

IS S-iuOD 

iSSDrioCD 


IADI = 2 IJ- 
lADIz? IJ: 
IADI=2 IJ- 
lADl-2 ijr 
lADIr? Ijz 
IAD1=2 
I ADi= 2 
IADI=2 
lADIrZ 
IADI-2 
I AdI- 2 
IADlr2 
lADI z2 
lADI-2 
lAOIz? 
lADIr? 
lADlzZ 
lADl-2 


4 ^ JSS — ^uCiJ 
TSS-lQnn 

iss=idd5 

jS£rlljQD 

ISSriGQO 

ISSzU'DG 
ISS-lLGC 
iss^icno 

ISS^lOOci 
TSS::lDaC 
iSS-lGoG 
ISS^IDOD 
ISS=1GDg 
ISS^IGCL 
j S S - 1 u 3 U 
ISSriGDD 
ISSulDOD 
ISS-I03D 


iSSG'iCCC 

ISSO-iCC? 

iSSDrlGCO 

1SSD=1dc2 

ISSDriOCiO 

ISSCrlCQD 

issr^iroc 

ISSd'IOCO 

ISSLriOGO 

ISSDriCOG 

] ssg=idoc 

IS£d=1CCO 

ISSDriGOC 

ISSDrlPGG 

ISSO^IOQO 

ISSDzicCD 

ISSDrlOOO 

ISSLrlGDG 


i 

5 

6 

iJz 7 
1J= 8 
IJ- 9 
IJzlu 

ij-n 

IJ-12 
I J-13 

1J = 15 
I J = 16 
lJ = l7 
iJrie 
lJ-19 



00 

00 


HLIRIC 


1 

X 

Y 

FT. (NHl ) 

FG( NH2) 

FGCNH!) 

?0 

9.432-qI 

1. 03C-*30 

9 .020-01 

9.320-31 

1 .000*30 

19 

9.432-01 

9.035-01 

9.823-01 

9.320-31 

i .000*30 

18 

9.432-01 

557-01 

9.821t01 

9.3 21-tEU. __ 

..l.Q 00 3.30 

17 

9.432-01 

9.228-01 

9.622-Cl 

9.522-01 

1 .000*00 

16 

9.432-01 

8.821-01 

9,626-cl 

9.526-01 

1 .000*03 

IS 

9.432-01 

8. 321-CI 

9.832-01 

9.532-31 

1 .000*00 

lit 

9.432-01 

7.721-01 

9.042-01 

9.542-01 

1.000*30 

13 

9.432-01 

7.D27-Q1 

9.858-01 

9.553-01 

1.000*03 

12 

9.4 32-01 

6.253-01 

9.882-01 

9 .5 3 2.r3J 

J ..000*30 

11 

9.432-01 

5.424-01 

9.914-01 

9.914-01 

1.000*30 

10 

9.432-01 

4.576-01 

9.957-01 

9.957-01 

1 .000*30 

9 

9. 432-01 

3. 747-01 

1 .col* on 

1 .3 ol*oo 

1 .000*30 

8 

9.432-01 

2. 973-01 

l,007« 00 

1.307*00 

1,000*00 

7 

9.432-cl 

2. 279-01 

1.013^00 

1.313*30 

1 . 000*30 

6 

9.432-01 

1 .679-01 

1 .Ol?* OQ . 

. I._119320 

1 .000*00 

S 

9,432-01 

1. 179-01 

1 .025*00 

1 .325*30 

1 .000*00 

4 

9.432-01 

7. 716-02 

1.031fC0 

1.331*00 

1 .000*00 

3 

9.432-01 

4. 47B-U2 

1 .036* 00 

1.336*03 

1 .000*30 

2 

9.4 32-qI 

1.949-02 

1 .040*00 

1.340*00 

1 .000*00 

1 

9.432-31 

0. 003 

1.043*00 

1.343*00 

1 .000*30 


dat^ 


F3<Nh12) 

F0iNH2i 1 

FS iwHM » 

FGINH32) 

BETA 

6,5 55-04 

-3.021-02 

o.dod 

-0,000 

3.417-Ofc- 

1.649-03 

-3.820-02 

0.000 

0,000 

5.03S-02 

3..7 0 8-0 3. 

-3 .814-nZ 

0 . non 

n.DOQ 

1.1S6>DI 

6. 554-03 

-3,799-02 

0.000 

0,000 

1.991-01 

1.306-02 

-3 .769.02 

0.030 

. 0.000 

i*ni5=Di- 

1,446-02 

-3.714-02 

0.030 

0.000 

4.310-01 

1.991-02 

-3.619rQ2. 

. 0,000. , __ 

— 0*.00il 

S.BllrOL- 

2.655-02 

-3 ,465-02 

0. 000 

0,000 

T.518-01 

.3 . _4 58-02 

-3.225-02 

0.000 

0.000 

9.Ta7>01 

4.4 14-02 

-2 .869-02 

0*000 

0,000 

1.119*00 

.5, 536-02 

-2.364-02 0.030 

--D.o.ao 

1.29Q«oo 

6.827-02 

-1 .681-02 

0.000 

0,000 

1.435«db 

0.267-02 

-8 .q76-03_. 

_ja.j3.ao 

_Q..OOQ 

1.543*00 

9.0 1 3-02 

2 .384-03 

0.030 

0,000 

1.609*00 

J. l4Qrfll 

1.409-02 

n.nan 

Q.DOQ 

1.637*00 

1.294-01 

2 .635-02 

0.000 

0.000 

1.632*0o 

1.43 9-01 

3.84lr02 

o.*ooo 

_fl«-0(10 

1.606AQO 

-1.58 1-01 

4.983-02 

0. 000 

0.000 

1.566*00 

-1. 76 0-01 

6.20,7-02- 

Q.oaa 

Q*oao_ 

ItJLMiOa- 

-1.4 19-01 

7.517-02 

D. O 30 

0,000 

1.302*00 



JX= 2, VARIABLE AT X( 2 )= 





level 

1 

l)-vel 

^ 4 4c 4c 



IZr 

1 

2 

3 

4 

5 

6 


Z - 

.3000 

» 1 0 C u + 0 0 

. 23Q0 + D0 

.3QD0+0C 

-«-4-30n+Di3 

-,.^011+013 

20 

.1300+01 

.1030+01 

• 1 L03+ 01 

.1330+31 

. lDDO+01 

.1300+01 

.1000+02 

19 

.9B10+C0 

.1000+01 

.1003+01 

,1330+Dl 

.1300+31 , 

. .laoo+ai — 

. 1000+01 

IB 

.9561+00 

.1030+01 

. 1C03+01 

.1330*31 

. 1300+01 

.1300+01 

. 1000 + 01 

17 

.9243 + 00 , 


..•.10QD..+ 01 .. 

. .1300 + 31,, 

..1DD.D+DJ^, _ 

..aiHQn+ni 

. lono+Di 

16 

.8842+00 

.1000+01 

• 10 03+01 

.1300+31 

. 13 DO +31 

.1300+01 

• lOuO + 01 

15 

.8349+00 

.luCD+Dl 

.1030+01 

.1330+01 

.1000+01 . 

. 13 0_D+_CL1 

_i..1Dd.O+01 

14 

•7759+GC 

.9999+30 

.9999+30 

.9599+03 

. 99 95 +30 

,9999+00 

.9999+00 

13 

.7375+00 

.9959+00 

.9999+00 

.9999+30 

. 9999 +30 

. 5 9 9 9 * n n 

.9999+00 

12 

.6310+go 

.9999+00 

.9999+00 

.9999+00 

.9999+00 

.9999+00 

. 9999+00 

11 

.5489+00 

.9998+00 

. «599B + 00 

,999R+nn 

t 95.98 +D .0 . 

,9998+00 

•9998+00 

10 

.4645+03 

.9997+00 

.9997+00 

.9997+00 

.9997+00 

.9997+00 

, 9997 +QO 

9 


.9997+00 

.9597+00 

.5997+00 

.9997+30 

.9997+00 

.9997*00 

8 

.3040+00 

.9996+00 

,9995+00 

•^ 996+00 

.9996+00 

.9996+00 

. 9996*00 

7 

.2339+qq 

.9996+00 

.9 996+ CC 

..9? 9 6 + 30 

. 99 9 6+00 

,9996 + nn 

t 9996*00 

& 

.1730+00 

.9996+00 

.9956+00 

.9996+00 

. 9596+00 

,9996+00 

.9996+00 

5 

.1218+00 

.9956+00 

.9955 i03 

.9996+00 

. 99 96 + go 

. 9995+00 

. QQ96+00 

4 

.7993-01 

.9996+30 

. 9 9 9 5 + 0 r- 

.9996+00 

.9996+00 

. 99964 OD 

.9996+00 

3 

.4653-02 

,9995+00 

.9596+00 

,9996+00 

,999b*_0D .. 

.9_9,9_6+.DI1 

.9996+00 

2 

.2327-01 

.1001+01 

.1001+01 

.1331+01 

. 1001+01 

.1001+01 

. 1001 + 01 

1 

.0000 

.7161+00 

. 7 1 e 1 + 03 

.7181+00 

. 7181 +0D 

.. 7.L81 + D_Q_ 

-^-7 18 1 + 00 


00 

VO 



VO 

o 


JX- 2t VARIABLE Xt i)~ .9H7 21 + CG 





LFVE L 

1 i?! ♦♦ 

V-VEL 




17- 
Z = 

1 

.ODDO 

a 

. 1C:GD + 00 

3 

.2300+00 

4 

• -CnG-^DD 

5 

.4000+00 

6 

. 5000+00 

20 

• 1 D D D U 1 

-.4918 -q4 

-.4916-04 

-.4916-04 

-.4916-04 

-.4918-04 

- • 4918-04 

19 

.98 IG+Ou 

- .3791 -D3 

-.8791-03 

-.8791-C3 

-.8791-03 

-. 8791-03 

-. 8791-03 

18 

• 9 5 6 1 ♦ L' G 

-.2318-02 

- . 2016-Q2 

-.2018-32 

-.2018-32 

-.2318-02 

-.2018-02 

17 

. 9 2 4 3 + D C 

-.3473-02 

-. 3473-02 

-.3473-02 

-.3473-32 

-.3473-02 

-.3473-02 

16 

. B8M2 + GC 

- .5296-02 

- .5296-02 

-.5296-Q2 

-.529^-02 

-.5296-02 

-.5296-02 

15 

•B349+CG 

- .7520-02 

-.7520-02 

-.7520-02 

-. 7520-02 

-.7520-02 

-.7520-02 

14 

• 7 7 5 5 + D Cl 

-.1D15--1 

IQIS-Cl 

-.1315-31 

-.1015-31 

-.1315-01 

1015-01 

13 

.7o754QC 

-.1312-01 

-.1312-01 

-.1312-01 

-. 1312-01 

-.1312-01 

-.1312-01 

12 

.631D“»-Du 

- .1631 -01 

1 6 31 - 

_..“.1^31''01_ 

-.1631-01 

-.1631-01 

-. 1631-01 

11 

.5 489+ GO 

-.1952-01 

1 952-01 

-.1952-01 

-. 1952-31 

-.1952-01 

-. 1952-01 

10 

-4645+00 

-.2261-01 

-.2251 -01 

-.2251-01 

-. 2251 -01 

-.2251-01 

-.2251-01 

9 

•381B+LQ 

-.2-33-Ci 

-.25X3-01 

-.2503-01 

-.2503-01 

-.25o3-ol 

-.2503-01 

8 

i3040+G0 

-.2691-01 

-.2691-01 

-.2591-01 

-.2691-01 

-.2691-01 

-. 2691-01 

7 

•2339+GD 

-.2805-01 

-.28c5*'Dl 

-.23 05-01 

-.2805-01 

-.2805-01 

-.28o5-D1 

6 

,1730+OG 

-.2852-01 

-.2852-01 

-.2352-0 1 

-.2852-qI 

-.2852-01 

-.2852-01 

5 

•1216+UU 

-.2842-01 

-.2642-01 

-.2“> 42-01 

-. 2842-01 

-.2842-01 

-.2642-01 

4 

.7993-01 

- .2797-01 

-.2797-01 

-.2797-01 

-.2797-01 

-.2797-01 

-.2797-01 

3 

• 4650-01 

-.2670-01 

-.2673-01 

-.2570-0 1 

-. 2670-0 1 

-.2570-01 

-. 2670-01 

2 

.2027-01 

-.2589-01 

-.2589-01 

-.2589-0 1 

-.2589-01 

-.2 5 89-ill 

-. 2589-01 

1 

• 0 3 D 0 

.3 DOG 

• n c 0 0 


.0000 

.3300 

.0000 





JX= 2f variable: at X( 2>r .9M321 + 0D 





LEVEL 

1 

W-VEl 




lZ~ 

1 

? 

3 

4 

5 

6 


Z r 

.0000 

. IGOu+uD 

.2^00+30 

. 30DO+QC 

.430C+D0 

.5DDD+Q0 

20 

• 1 0 C 0 + Pi T 

. 0 0 0 Q 

. 2C1 &-GS 

.1346-38 

-.36 16-38 

- f 5 5 4 .lr_H 8 _ 

-.6183-08 

19 

.98lD+ca 

.0000 

. 1 397-oa 

,97 28-39 

-.2746-08 

-.4653-08 

-.5588-08 

18 

•9561+00 

•^000 

-.4657-09 

-.1475-39 

-.1359-09 

-.l?9P-o8. 

. _:i^4l9l-D8 

l7 

.9243+00 

.ocao 

-.2375-08 

.6750-39 

-.2995 -08 

-.37 19-09 

. 1397-08 

16 

.3842+00 

.0000 

.79ifa-08 

-.2569-38 

liP?rD 8 

.1386-08 

. 1863-08 

15 

.3349+00 

.0000 

.4191-06 

.7366-36 

-.2434 -08 

-.1301-08 

-.4657-09 

14 

.7759+00 

.OGDG 

.4191-38 

.2428-38 

-. 3775-08 

-.3095-08 

-.2794-08 

13 

.7 075 + 00 

.0030 

.5122-08 

-.2330-38 

-.39 10-09 

.3529-09 

. 1397-08 

12 

.5310+00 

.0300 

.0003 

-.-'5Q2-3& 

. e4 43-09 

.54 12-p8_. 

.6985-08 

11 

.5489+00 

• 3 0 0 

.8848-08 

-.9947-39 

-. 1722-38 

.5486-08 

.6985-08 


,4645+00 

.ODorj 

-.5586-08 

-.4476-38 

18??-a8 

.7322-09 

.9313-09 

9 

.3818+00 

.3000 

-. 1 397-C6 

-.1363-38 

. 4S73-39 

-.2653-08 

-.2328-08 

8 

• 3Q43 + LiQ 

.3330 

.2794-08 

.4980-39 

-.2063-08 

*18 4 7-Q9._ 



7 

. 2 33«, + Q0 

• 3030 

-. 1164-C7 

.3333-n6 

. 5206-08 

-.5636-08 

-.9779-08 

6 

.1730+00 

. 3 G J 0 

-.5iZ2-Gr 

-.9750-39 

-. 1453-08 

.3 04 2-08.. 

.3725-08 

5 

• 1218 + (]Q 

.3000 

-.3725-05 

.4593-38 

. 2582-08 

-.1252-08 

-.1863-08 

4 

.7993-01 

.3030 

.4657-09 

.2311-36 

-.7713-08. 

-.5051-08 

-.6985-08 

3 

.4650-Cl 

,3030 

-.9313-09 

.2736-08 

. 2355-38 

.3589-08 

.5588-08 

2 

.2u27-D1 

. 0030 

.5588-08 

.1377-39 

- , 54 65-08 

?5 9 2_-Q9 , 

.46AIr09_ 

1 

.DDCO 

.3030 

- . 9 3 1 3 - fj9 

.2736-38 

.2355-08 

.3589-08 

. 4DDD-q8 



VO 

to 


JX = 

2, VARIABLE 

AT X( 2>= 

.94321+00 








level 1 


RhO 




IZ = 

1 

2 

3 

4 

5 

6 


Z - 

.0000 

. ICCQ+QQ 

.2030+00 

• ^OCG+OD 

.4300+00 

•5Q00+00 

2D 

•1 QOG+01 

.1000+01 

.1000+01 

.1000+01 

. lODQ+01 

.1300+01 

. lOOD+01 

19 

•9810+00 

•1 CQO+Ol 

. ICOO+Cl 

.1000+01 

. IQDD+Dl 

. . ..IIQ.O.+JIL 

. lODQ + Dl 

18 

.9561+00 

.1000+01 

. 1G00+ 01 

.1000+31 

. IDGD+Ol 

.1000+01 

. 1000401 

17 

•9243+uU 

.1000+01 

.1003+01 

. 1000+01 

. 1000 + 01 

.1000+01 

. 1000 + 01 

lb 

•3342+uD 

.1000+01 

• 1 QCO+ 01 

. 1000+01 

. 1000+01 

.1000+01 

• lOOO+Ol 

15 

.B349+Q0 

.1000+01 

. 1000+02 

.1000+01 

. iDoo+ai . 

.1000+01 

. 1000 + 01 

14 

.7759+QG 

.1000+01 

. lODO+Ql 

.1000+31 

. lOOD+Ol 

,1000+01 

. 1000 + 01 

13 

.7075+00 

•1 coo+oi 

. 1 DDO+ Q1 

.1000+01 

. 10 00+0 1 

.ipoo+o_i 

. 1000 + 01 

12 

.6310+00 

.1000+01 

. 1000+ 01 

, iODD+01 

. 1300 + 01 

.1000+01 

•1000+01 

11 

.5469+00 

.1000+31 

.1000+01 

.1000+01 

. 1000 + 01 

.1300+01 

■ IDOO+Dl 


.4645+00 

. 1 QDO+01 

.ICDO+Dl 

. 1000+01 

. 1000 + 01 

.1000+01 

. 1000 + 01 

9 

.3818+00 

.1000+01 

. 1 0 D 0 + 0 1 

.1000+01 

. lono+oi 

.1000+01 

. lOOn+ni 

8 

.304 o+GO 

.1000+01 

. lDDO+01 

. lOOD+31 

. lODD+01 

.1000+01 

. lDDO+01 

7 

.2339+00 

.1000+01 

.1000+01 

. 1000+01 

. 13DD + 01 

.1000+01 

. 1000+01 

6 

.1730+00 

•1 000+01 

. 1 0D0+ 01 

.1030+01 

. 1000+01 

.100 0+01 

. 1000 + 01 

5 

.1216+00 

.1000+01 

. lGOO + 01 

.1000 + 0 1 

.JOILCiJ+M. 

.1000+01 

. 1000+01 

4 

.7993-01 

.1000+01 

. 1Q0Q+ Cl 

. 1000 + 01 

. lQOO+01 

•1300+01 

• iQOO^Ql 

3 

.4650-01 

.1000+31 

.1000+01 

. 1000 + 0 1 

. 1000+31 

.1300+01 

. IDDQ + Dl 

2 

.2027-01 

.1023+31 

. 1 023+ G1 

. 1023+01 

. 1023+01 

. 1 02 3 + 01 

. 1025+01 

1 

.Doro 

.1339-02 

• 1 309-02 

.1309-02 

. 13 09-02 

.1310-02 

.1310-02 



JX- 

2, VARIABLE 

aT X( 

2)- .94321+00 








LEVEL 1 

(t 

ENERGY 

♦ . 

— 


12- 

. 1 

•n 

C 

3 

4 

1 

1 

1 

i 

! 

6 


Z - 

.3030 

.1003+30 

.230D+3C 

... • ?aDa-+0G-.-. 

.43GD+DD_ 

. SDCD^DQ 

20 

.1300+01 

.3030 

.0003 

.0330 

.0300 

.3300 

• 0000 

19 

.931u+D0 

.0030 

. Li 0 J 

.0300 

.03 00 

..31QQ 

... _ 

18 

.9561 +QD 

.0030 

. 0 L C 3 

.0330 

.0300 

.3300 

. 0300 

l7 

.9245+00 

. 3330 

.3000 

. 0303 

• 3000 

.. .3:200.- 

. ....OQilIL 

16 

.3(342 + 00 

.3000 

. n G n 0 

.3300 

.GOOD 

.3300 

• GOOD 

15 

.B349+LQ 

.3030 

.3c-00 

.0300 

.DODD 

.3300 

, GOOD 

14 

.7759+CD 

.3000 

. ODCJ 

.0300 

.0300 

.3300 

.0030 

13 

■ 73 7S + 0C 

.3000 

. n0D3 

. 0300 

• 0300 

.33 00 

..G.DQ.0 

12 

.6310+00 

.3030 

.3033 

.0330 

. 0300 

.3300 

. GO Q g 

11 

.5489+00 

.3000 

.0003 

. 0333 

. GOOD 

• 33QD 

,.0000 

10 

.4645+00 

.3000 

.0030 

.0300 

.0300 

.3300 

.0000 

9 

.361B+0Q 

•3000 

.0003 

.0330 

. 03 00 

. .Ojc-a 

.0000 

8 

.3043+30 

.3030 

.3003 

. 3303 

• 0330 

.3300 

.0000 

7 

.2339 +UD 

. 3 0 J 0 

.0003 

.0303 

. nono 

.3300 

. . -.-•-DU D.D 

6 

.1730+00 

.3000 

.0030 

.0300 

. GOOD 

.3300 

.0000 

5 

.1218+00 

.3000 

• 0003 

.D3 03 

. OOGu 

.3300, 

.. ...nono- 

4 

.7995-01 

. 30 30 

. 0 u U D 

.0330 

. 0300 

.3300 

.0000 

3 

.4650-01 

.3330 

. 0 003 

.0330 

. ODGJ 

. 3.3 3 D 

. 0 0 D n 

2 

.2327-01 

.3000 

. 0 [. u 0 

. 3300 

.0300 

.3300 

.0000 

1 

.3000 

.3030 

. 0030 

.3330 

. D30Q 

.3300 

.. ...QQ.OQ 


VO 

L3 



VO 

-IN 


JX- 2, VARIABLE aT X( 2»- .94321+00 

LEVEL 1 



IZ- 

1 

2 


2 = 

.3030 

.Too 3 + 06 ' 

20 

.13GQ+U1 

.3122+00 

.3122+ 00 

19 

.9813+00 

.3122+00 

.3122+00 

18 

.9561+00 

.311 1+30 

.3111+03 

17 

.9243+CQ 

. 3 09 3 + 00 

• 3C93+ DO 

16 

,38q2+GD 

.3063+00 

• 3 Q 6 3 + 0 G 

15 

.9349+OD 

.3013+00 

. 30 13 TTT' 

14 

.7759+00 

•2936+QO 

.2936+ GO 

13 

. 7 O 75 + DO 

. 2 B 2 D + QV. 

. 2323+00 

12 

•6310+00 

.2648 + q(^ 

. 2 64 S + UQ 

11 

.5489+00 

. 2 4 Q 6 + 0 0 

. 2406 + 00! 

ID 

.4645+03 

.2074+30 

.2074+00 

9 

. 38 18+ QQ 

, 1640+00 

. 1 643 + 00 

8 

.3040+00 

.113 1 +30 

. IICI + CO 

7 

.2339+00 

.4718-01 

.4716-01 

6 

.1730+00 

.2223-al 

.2223-01 

5 

.1218+00 

. 9468-01 

.9468-01 

4 

.7993-01 

.1693+00 

.1693+00 

3 

•465u~Dl 

,2879+30 

.2679+00 

2 

.2027-01 

.6445+00 

.6445+30 

1 

. 3 3 0 j 

-.1178+01,. 

-.1178+ D_i, 






• 3qc!0*00 
. 3122+30 
. 3122+00 
. 3111+30 
. 3Q93+0C 
. 3 n63 + go 
.3013+36 
._2?36+,QD_ 
. 29 20 + 30 

• ?& 4 8 +p_p._ 
.24 o6+qq 
. 2D74 +3D 
."15 4 0+3 0 
. UD1 + 3D 
.4713-01 
.2223-01 

.-.?4 6_8_-Q„L 
.1593+33 

_.2ai9.+J3Q 
. 64 4 5 +0 0 


.4300+00 

.3122+00„ 

.3122+00 

• 3 1 il 10 0_ 

. 3093+30 

.3363+00 

.3313+00 

.2936+00 

.2B20+D0 

.2548 +00 

.2406+00 

.2074»DD 

.1640+00 

. 11 01+00 

.4718-01 

.2223-01 

.9468-01 

.1693+00 


.5000+00 
. 3122+00 
.3122+00 
. 3111+00 
. 3093+00 
. 3D63»D0 
. 3013 + 0(3 


. 2820+00 
.2648+00 
.2406+00 
.2074+00 
. 1640 + 00 
. 1101+00 
,47i8-0i 
.2223-01 
.9468-01 
• 1693 + 00 
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.6717-04 
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.6717-04 
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.5716-04 
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.6716-04 
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.6716-04 
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.6716-34 

.5716-04 
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13 

.7075+00 

.6715-04 
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.6715-34 

.6715-34 

.5715-04 

.6715-04 

12 

.6310+00 

.6715-34 
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.6715-04 

.5715-04 
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11 

.5489+00 
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.6713-04 

.6713-34 
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.5713-04 

.6713-04 

10 

.4645+00 

.6712-04 
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9 

.3818+00 
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.3340+00 
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19 
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.9561+00 
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.9243+00 
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.8842+00 
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.3349+00 
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. lOGO+O 1 
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• lOOg+01 

11 

.5489+00 

.luQO+Ql 

. 1003+01 

.1300+01 

. IDGO+Dl 

.1000+01 

. lODD+01 

10 

.4645+00 

.1 CD 0 + 01 

. 1CC0+ 01 

.1300+01 

. 1300 + 01 

.1000+01 

. IDDD + Dl 

9 

.3816+00 

• 1 CD 0 + 0 1 

. ILCO+Cl 

.1300+01 

. lODO+Dl 



. IQDD + Ql 

6 

.304D+DD 

.lODD+01 

. ILCD+Cl 

.13D0+gl 

• IQOC+Dl 

.1000+01 

. 1000 + 01 

7 

.2339+00 

.1000+01 

. 1GDD+ 01 

.1300+01 

. loco+ai 

.1000+01 

.1000+01 

6 

.1730400 

. 1 DOD + g i 

. ICGD+Ul 

.1300+31 

. lOCD+Dl 

.1000+01 

. ICDD + Ol 

5 

.1218+00 

.1000+01 

.ICCD+Cl 

.1330+01 

. lDDD+01 

.1000+01 

. IDDO + Dl 

4 

.7993-01 

.1000+01 

• 1 003+ Oi 

.1300+01 

. 1000+01 

.1000+01 

• IDOD + Ql 

3 

.4650-01 

.1000+01 

.1003+01 

. 1300+01 

. IDOO + gl 

...r^Q0_0+01 

. 1000 + 01 

2 

.2027-01 

.9945+00 

, 9945 + DC 

,9945+00 

.9945+00 

.9945+00 

.9945+00 

1 

• 0 0 GG 

.1 872+qI 

.1872+01 

.1972+01 

. 18 72+01 

.1872+01 
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lADI-i 

IJ- Z 

iss-inoD 

ISSC^IDOO 

IADI-1 

IJ= 3 

ISS^IDOD 

ISSD^IDCO 

lADI^l 

iJ- 4 

ISS::lODa 

ISSD-IOOC 

lADIrl 

IJ- 5 

issncoo 

lSSD-1 DUD 


lADIr? 

IJz 2 

ISSr iQOQ 

ISSD-IOOO 

1ADI=2 

IJ- 3 

ISSzluDD 

iSSDrlOOa 

I AoI- 2 
IADIr2 

1J= 

I J- 5 

ISS - 1 [JQC _ 

IS S:: 100 0 

I SSC!?AD^0Q 

ISS 0=1 DOG 

lADIrZ 

1J= b 

ISSrlOnn 

ISSD=1D00 

IADI=2 

IJ= 7 

ISS-IUDD 

iSSDrlOOO 

lADirZ 

IJ= 8 

tSS-IOOD 

lsSD-1 nno 

lADlrZ 

IJ- 9 

ISS- lOOQ 

lSSD-1000 

lADI^Z 

I JrlO 

iSSrlCDD 

ISSDrlDOn 

IADI-2 

I J-li 

ISs=^1GlIQ 

ISSD-IDOO 

IADIr2 

I J-12 

ISS::iDQQ 

ISSd=10D0 

lAOIzZ 

I J = 13 

ISS-IuQq 

ISSDri DOQ 

IA0I=2 

IJ-14 

TSS-IDOO 

ISSDziDOO 

IADI=2 

I J-15 

ISS-IDOO 

ISSD-IOOO 

lADI-2 

I Jzlb 

ISS^IOOG 

ISSD-lOOn 

IADI-2 

1J = 17 

ISs-lOOD 

ISSO-IQOC 

lADIrz 

I J=ie 

ISS^ICOD 

issd=idqo 

IADI-2 

I J-19 

I5S-1UD0 

ISSD-IOOO 
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iETPlC DAU 


1 

X 

r 

rrHNHl 1 

NH2) 

FG< NH3) 

F3 (NH12» 

FG1NH21 ) 

F3(NH31) 

FCINH32I 

BETA 

20 

9.602-01 

1. 000*00 

9.8Q9-Q1 

9.309-01 

1.000*30 

-6.684-04 

-4.014-02 

0.000 

-- 0.000- - - 

... 3.17Qt06 


9. 682-01 

9. 805-01 

9.809-01 

9.309-01 

1 .000*03 

-1.677-03 

-4 .013-02 

0.030 

0.000 

5.277-02 

IB 

9.682-01 

9.552-01 

9 ,6 in-ni 

. 9.3lDr21 . 

. -1 .3Qfl.+aa... 

-3. 451-03 

-4 .Q07-Q2 

n* non 

n.ODD 

1 . 212 -ni 

17 

9.682-01 

9.228-01 

9.bil-01 

9.S11-D1 

1 .000*30 

-6.667-03 

-3.993-02 

□ . 000 

0.000 

2.087-01 

16 

9*682-01 

8.621-01 

9.815-01 

9,315-01 

1 .000*00 

-1,324-02 

-3.966-02 

D.OoO 

O.OQQ . 

-_..3.ia3-Ql 

15 

9.682-01 

8.321 -01 

9.321-01 

9. 821-01 

1 .300*30 

-1.471-Q2 

-3.915-02 

0.000 

0.000 

4.522.01 

m 

9.682-01 

7.721-01 

9.832-01 

9,332.31 

1 .300*30 

-2.327-02 

-3.B2a-:Q2 

0,030. 

- G*ooa — 

6..-lQ6 = Qi. 

13 

9.682-01 

7. 027-01 

9.e«0-ol 

9.348-01 

i .000 ,00 

-2.706-32 

* 3 « 6 6 6 — fi 7 

0. 000 

O.OOQ 

7.906-01 

12 

9.6^2-iu 

. 8i2SirQl 

9.^72-01 9.372-31 

1 .000*30 -1.5Z9-02 

-3.464-02 

Q.QQQ 

O.QQO 

9^85 j-oj 

11 

9. 682-01 

5.424-01 

9.905-01 

9.905-01 

1 .000*00 

-4.513-02 

-3,135-02 

0.000 

C.OOO 

1^183*00 

Id 

9.682-01 

4. 576-01 

9,948-01 

9. 940. 01 

1 .000*00 

-5.625tQ2 

t2. 667-02. . 

. O.ODD - 

-o*ooa — 

.-U.32HQ0 

9 

9.682-01 

3. 747-01 

1. 000*00 

1.303*00 

1 .000.00 

-7.321-02 

-2 .03i-02 

□ .ODD 

O.OQQ 

1.534*00 

8 

9.682-01 

2.973-01 

1*006*00 

1.306*30 

1 f0O0*O0 

-6.532-02 

-1 .2U-Q2 

Q.QOD . 

.Q.QQQ 

1 . 6621 QQ 

7 

9.682-01 

2.279-01 

1 . 013* QQ 

1.313*30 

1 .000*00 

-1.317-01 

-2.369-03 

□ •030 

O.OQQ 

1.749*00 

6 

9.682-01 

1.679-01 

1 . l19«00 

1 .319.00 

1.300*30 

-1 . 1 8 5-n 1 

6 . ,.30-03 

0.030 

O.OQQ 

1.798*00 

5 

9.682-01 

1.179-01 

1 .025*00 

1.325*00 

1 .000*00 

-1.351-01 

2.020-02 

0.030 

O.OOQ 

1.815*00 

9 

9.682-01 

7.716-02 

1.031* 00 

1.331*00 

1 ,000*00 

-1,507-01 

- 3,162to2 

Ot-ODQ _ 

- Q.QQQ 

__J.8Q6.tDQ 

3 

9.682-01 

4.478-02 

1 .036.00 

1.336*00 

l.ODD^OO 

-1.671-01 

4.254-02 

0.000 

0.000 

1.785*00 

2 

9*682-01 

1,949-02 

1 • 041 ♦ 00 

1,341*00 

. l,3QD*Qa_ 

t1. 971-01 . 


. 0*000-- - - 

-I1.QQQ 

_.1.7Jl*QQ. 

1 

9*682. q1 

0.030 

1.044*00 

1 .344*00 

1 .000*00 

-1.693-01 

7.060-02 

0*000 

0.000 

1.51B*nn 



o 

4N 


JX = 

3t VARI/SSLE 

ftT X ( 3 jr 

.96821+00 












LEVEL 1 


U-VEL 




12 - 

1 

2 

3 

4 

5 

6 


1 - 

.0000 

. 1 JD3 + QC 

. 23_0_0 +D 0 

. 3rmn+nn 

.4300+00 

*5000+00 

20 

.9995+00 

.1300+01 

• 1 C03+ 01 

.1300+01 

. 1000+01 

.1300+01 

• 1000+01 

19 

• 95 04 + [;n 

. 1000+31 

. 1300+01 

. 1300+31 . 

. 1000+01 

-.1 0 00 + Pi 

. 1 nnn+ni 

le 

•9555+00 

•lOQO+jl 

. lODO+Ql 

. 1300*01 

. lDOO+31 

.1300+01 

. 1000 + 01 

17 

.9237+00 

. 1000+01 

. 1 GC 0 + 01 

.1300+01 

. lOGO+Ql 

. 1000+01 

. 1000+01 

16 

.8837+00 

. 1000+01 

. 1000+01 

.1300+01 

. 1000+01 

• 1 000+01 

. 1000+01 

15 

.8345+00 

.1030+01 

. IDD3+D1 

. l3DDtO_L__. 

-■tlQQ[I.+JII_- 

.1300+0.1.. . 

. 1000+01 

14 

.7756+00 

.9999+30 

, 9 9 o 9 + UD 

.9?99*00 

.99 99+00 

.9999+00 

,9999+00 

13 

. 70 72 + 00 

.9999+00 

.9999+00 

• 9? 99 + 0 n 

>9999 +00 

.9999+00 

.9999+00 

12 

.6307+00 

.9999+00 

.9999+00 

.9999+QQ 

.9999+00 

.9999+00 

.9999+00 

11 

.5487+CO 

.9996+Ou 

.9998+00 

. 9998+00 

.9998+00 

• 99 98+00 _ 

.9998+00 

ID 

. 4 645 + uC 

.9997+03 

•9997+00 

.9997+00 

.9997+00 

.9997+00 

.9997+00 

9 

.381&+J0 

.9995+03 

.9996+00 

.9996+00 

.9996+00 

.9996+00 

.9996+00 

b 

.3040+00 

.9996+03 

.9996+00 

,9996+00 

. 9995+00 

.9996+00 

.9996+00 

? 

.2339+00 

.9995+03 

.9995+00 

.9995+00 

. 9995+00 

.9995+On 

.999ci+D0 

6 

. 1 7 30 + 00 

•9995+00 

.9995+00 

.9995+30 

.9995+00 

.9995+00 

.9995*00 

c; 

.1219+00 

.9995+03 

.9995+00 

.99 g5 + 00 

. 9995+00 

. 9995+00 

•9995+nn 

4 

.BOOG-qI 

.9994+00 

.9994+00 

. 9994+00 

• 9994 +00 

. 9994+00 

•9994+00 

3 

.4655-01 

.9994+00 

, 9994+00 

.9994+00 

.9994+00 

.9994*00 

.9994+00 

2 

. 20 3u -G1 

.9977+00 

.9977+00 

.9977+30 

. 9977+00 

.9977+00 

.9977+00 

1 

• j j DO 

.7183+30 

.7183+00 

•Jl 8.3 4^0 

. 7183+00 

.7183+0n 

,71q3.+Q0 
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JX= 3, VARIABLE aJ X( 3)- .96821+00 





LEVEL 

1 

V-VEL 




IZ = 

1 

2 

3 

4 

5 

6 


Z = 

.0000 

. 1000+ 00 

.2000+00 

• 3000 + 00 

.4300+00 

.5000+00 

20 

•9995+00 

-.3579-04 

-.8579-34 

-.8576-04 

-.85 79-0 4 

-.8579-04 

-.8579-04 

19 

• 9304 +DD " 

— .9265-33 

-.9265-03 

-.9265-33 

-.9265-33 

-.9265-03 

-.9265-03 

18 

•9555+00 

-.2119-02 

-.2119-02 

-.2119-02 

-.21 19-02 

-.2119-02 

-. 2119-02 

17 

•9237+00 

-.3641-02 

3^-41-02 

-.3641-02 

-. 364 1-02 

-.3641-02 

-.3641-02 

16 

• B8 37 + 00 

-.5555-02 

-.5555-0? 

-.5555-02 

-. 5555 -02 

-.5555-02 

-.5555-02 

15 

.8345+00 

-.7892-02 

-.7692-02 

-.7392-02 

-.7892-02 

-.7392-02 

7892-02 

m 

•7756+00 

-.1D06-01 

- . 1 L 6 6 - 0 1 

-.1066-0 1 

1066-01 

-. 1066-01 

-. 1066-01 

13 

.7072+00 

-.1380-01 

-. 1380-01 

-. 1380-0 1 

1380-01 

-.1380-01 

-.1383-01 

12 

.6307+00 

-.1720-qI 

-.1720-01 

-.1720-01 

-.1720-01 

-.1720-01 

-.1720-01 

11 

.5487+00 

-.2065-01 

-.2 2.-6 5-01 

-.2065-01 

-.2065-31 

-.2365-01 

-.2065-01 

10 

.4645+00 

-.2593-01 

-. 2393-ul 

-.2393-01 

-.2393-01 

-.2393-01 

-.2393-01 

9 

.38 iB+oo 

-.2673-01 

-. 2b7B-Dl 

-.2673-01 

-.2678-qI 

-.2678-01 

-.2678-01 

8 

.3040+00 

-.2902-01 

-.2902-01 

-.2902-01 

-. 29 02-01 

-.2902-01 

-.2902-01 

7 

.2339+00 

-.3055-cl 

-.3 05 5 -01 

-.3055-01 

-.3055-31 

-.3055-01 

-.3Q55-D1 

6 

• l7 30 + 00 

-.3141-01 

-.3141-01 

-.3141-01 

-. 314 1-Ql 

-.3141-01 

-.3141-01 

5 

.1219+00 

-.3167-01 

-.3167-01 

-.3167-31 

-.3167-01 

-.3167-01 

-.3167-01 

9 

.8000-01 

-.3158-01 

-.3153-01 

-.3158-01 

-.3158-01 

-.3158-01 

-.3158-01 

3 

•4655—01 

-.2909-01 

-.2909-01 

-.29Q9-01 

-.2909-01 

-.2909-01 

-.2909-01 

2 

.2030-01 

-.2757-01 

-. 2 757-01 

-.2757-01 

-.2757-01 

-.2757-nl 

-.2757-01 

1 

• ooco 

.0030 

n /*' ^ r- 

• w 

.0100 

. OQDO 

.33 00 

• 0000 



o 

ON 


JX- 3, variable: at X< 3 >r .96821+00 



IZ- 

1 . 

LrVfL 

2 „_, _ 

1 

3 

W - vEL 
4 

5 

6 

— 


Z - 

.ODOC 

. 1000+ GO 

•z GOO * 00 

. 3000+00 

,4300+00 

.5000+00 


2D 

. 9995 + CD 

.0000 

• lt2b — u7 

. 72 5 9r.Q8. 

- . 1235.-J17._, 

-.1429-07 

1494-07 



19 

.9804+00 

.0000 

,7451-08 

,5?5 3-38 

-.9003-38 

-.111 5-07 

-. 1164-07 


18 

.9555+QU 

• 0 D Q 0 

6985-03 _ 

..._-..75 4 3-a9 _ 

. 1036-08 

-.17 19^08_ 

-. ^79^-a■8■■. 



17 

.9237+00 

.0300 

-.4657-08 

.5390-38 

.418^-09 

-.4313-08 

-.6985-08 


16 

.8837+00 

.0000 

. 1956~(i7 

-.TByO-Do 

-. 7050-08 

.2202-08 

.4657-08 


15 

.8345+00 

ODQQ 

. 1537-07 

.2522-07 

-. 1092-07 

-,i5&3-d7 

1863-07 


m 

•7756+00 

• OPOD 

. 1 72 3.-07 

„___!_539 5-a8 _ _ 

1661-07 

-.6370-08 

-.5122-08 


13 

.7072 + L(Q 

.0000 

.9779-08 

-.2373-08 

-. 1028-08 

.1457-08 

. 1863-08 


12 

• 6307 + 00 

.ODOD 

. 6519-08. 

_._“^§.5.55^0.8 

-.48 qD-08 

.1357-08 

.1397-08 


11 

.5487+00 

.0000 

. 2654-07 

-.1345-08 

-. 3974-08 

.7372-08 

. 1071-0? 


1C 

•4645+qg 

• 30 0 0 ..,.. 

. 7916-08 

-.6391-08 

-. 3437-08 

.7275-08 

.1211-07 


9 

.3818+00 

.0000 

-. 1071-07 

-,1i55-d8 

-. 1531-08 

-.9333-08 

-. 1211-07 


e 

.3040+00 

.3000 

-.4657-08 

-.6117-08 

3747-38 

-.4361-09 

,93] 5-d9 


7 

,2339+QQ 

. 0000 

-.2142-07 

.91 65-1 D 

156q-o9 

-.1441-07 

-.2142-07 


6 

.1730+03 

.0000 

-.1583^07. 

.. . - *.8 9 2 l-Jl 8 _ 

. . . 32iLl.^A. . 

.1 031-Ofl 

.4657-09 


5 

.1219+00 

• 3Q0Q 

-. 1909-07 

-.1500-06 

.8104-08 

-.256[j-d8 

-.7916-08 


4 

.BOQO'Ol 

.0000 

.7451-08 

-.1958-08 

-.1390-07 

-.1004-07 

-.8848-08 


3 

.4655-01 

.ObOD 

, 16i6-g7 

.1385-07 

.2653-08 

.3576-08 

•326q-o8 


2 

•203o-al 

.3000 

.2049-07 

.1406-07 

-.2309-08 

-.6971-08 

-.‘>313-08 


1 

• 3 0 0 0 

• 0 030 

.1659-07 

.1381-07 

.4187-08 

.3845-08 

. 3731-08 
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JX- Zf VARIABLL AT XI 3)= ,96821+CC 

level 1 



IZ- 

1 



Z r 

.0000 

. icno + DC 

20 

.9995+00 

,100G+Gi 

. 1 CO 3+ 01 

19 

.9804+Uu 

.1000+01 

.1003+01 

16 

.9555+00 

.9999+30 

.9999+00 

17 

.9237+00 

.9999+00 

.9999+00 

lb 

.8837+00 

,9999 + OC. 

. 9999 + 00 

15 

,8345+00 

.9999+00 

,9999+00 

14 

.7756+00 

,9999+OU 

. 9 9Q9+ 00 

13 

.7072+00 

.9999+00 

. 9 99 9 + no 

12 

,6507+u0 

.9999+30 

.9999+00 

11 

.5467+00 

.9999+00 

,9999+00 

10 

.4645+00 

.9999+00 

,9999+ 00 

9 

.3818+00 

.9999+DC 

.9999+00 

6 

,3040+QQ 

.9998+00 

.9998+03 

7 

.2339+00 

,9998+00 

.9998+00 

6 

.1730+00 

,9997+00 

.9997+03 

5 

.1219+00 

.9997+00 

.9997+00 

4 

.3000-01 

.1031+01 

,1001+01 

3 

.4855-01 

.1331+31 

. 1001 + 01 

2 

• ? 3 3 0 ■ 0 1 

.1042+31 

. 1042+01 

1 

• ODOQ 

.5721+30 

• 5721 + 00 






.23Q0+0D 

,1300+01 

. 130 0*01 

,9999 + [)3 

• 99 99 + Dp 

,9999+30 

,9999+00 

,9999+00 

.9999+00 

.9999+00 

.99 9 9 + PP_ 

.99 9 9+ on" 

,9999f00 

.9993+00 

. 999B+0C 

,9997+00 

.9997+00 

.1301+31 

.1301+01 

. n^^+oi 

.5721+00 


• 3000 + 00 , 
. lOQO + 31 

• 10 00 +01 
,9999+00 
. 99 99+00 
,9999+OQ 
. 9999+00 
.9999+30' 
. 99 99 +0i^ 
.9999+00 
. 999 9+0 0 
.9999+00 
.9999+00 
. 9998 +00 
. 9998+00 
. 9997+30 

• lOOl +0l 
. 10 01+01 

• 1042+31 
.5721 +00 


• Aooc + ga _ 

.1300+01 

J .1 2 QD 1 Q 1 . 
.9999+00 
.9999+00 

• 9 999+00 
.9?99 + OjO_ 
.9999+00 
. 9999+DU 
.9999+DO' 
. 9999+0 0 
,9999+00 
,9999+ no 
.9998+00 
•??9_8 + qo 
.9997+00 
.9 997 +0 0 
.1301+01 

• 1301 .+ 01 _ 
. 1042+01 
.5721+00 


. 5D00+QD 
. 1000+01 
. IDDO + Dl 
, 9999+00 
. 9999 + PO 
. 9999+00 
.9999+00 
.9999+00 
. 9999+00 
,9999+00 
. 9999+00 
• 9999 + 00 
. 9 999+00 
.9998+00 
. 9998+00 
.9997+00 
.9997+0 0 
. lODl+Dl 
_* loxU - iiii - 

. 1042+01 
.5721+00 



H 

O 

CX5 


JX- 

3, Variable 

at X ( 3 ) = 

.96821+00 






- 


LEVEL 1 

***** 

ENERGY 




iZz 

1 


3 

4 

5 

6 


Z - 

•3000 . . . 

. lODO+DO 

.2000+00 

.3000+30 

/4 onn+nn 

. 4000+00 

20 

•9995+00 

.0000 

. 0003 

• GOOD 

.0000 

.3300 

.0000 

19 

•98q4+00 

.3000 

. 0000 

.00 30 

. 0000 

.3000 

.0000 

18 

•9555+OQ 

.0000 

.0000 

.GOOD 

. 0000 

.3300 

.0000 

17 

•9237+00 

•GDDO 

.0000 

.GOOD 



- «30qD 

-^Doni] 

16 

•3637+00 

.0000 

.0000 

.0030 

• 0000 

.0300 

.0000 

15 

•3345+00 

.3000 

•.00 DO... .. 

.0030 

. 0000 

.0300 

• ODD 13 

14 

•7756+00 

.3000 

.0003 

.0030 

. 0000 

.3300 

• GOOD 

13 

•7072+00 

.0000 

• 0 000 

.0003 . 

•00^ 

.3300 

.0000 

12 

•6307+00 

.3000 

. 0 C 0 3 

.33 00 

. GOOD 

.3300 

.0000 

11 

.5487+00 

. 0000 

.0003 

.0000 

.0000 

■0300 

■■ .0000 

IQ 

.4645+00 

.0000 

.OuOO 

.0000 

.0000 

.0300 

.0000 

9 

.3613+00 

• 0000 

.OCDU 

.0330 

. GOOD 

.Onnn 

.0000 

8 

.3040+00 

.0000 

. 0300 

• 0300 

. 0000 

.3300 

.0000 

7 

,2339+00 

. 0000 

. 0003 

.0300 

. 0000 

.3300 

.0000 

b 

.1 73D+QQ 

.0000 

.0003 

.GOOD 

.0000 

.3300 

• 0000 

5 

.1219+00 

• ODOD 

. DCQQ 

.3300 

. 0000 

.3300 

.0000 

4 

•eoco -01 

.3000 

.0003 

•0000 

.0000 

.3300 

.0000 

3 

•4655-01 

.3000 

.0003 

.3000 

. 0000 

• 3300 

• 0000 

2 

.2030-01 

.3000 

. 0 033 

.0000 

.0000 

.3300 

.0000 

1 

.0300 

.3000 

.0003 

.GOOD . 

.•_Q.Q_Q..o 

. Tinn 

.. + 0D.0.D 


W' 


JX= 3, VaRIaBLl at X( 3)- .9b621+D_0_ 





LEVEIL 

1 dt*** 

MUT 




12 = 

i 

c. 

3 

4 

5 

6 


Z - 

.3000 

. 1003+ GO 

.2300+30 

. 3Qoc+on 

.4300+00 

. 5DQO+00 

20 

.9995+00 

.3203+qo 

. 32 33*00 

.3203+30 

. 3203+30 

.3203*00 

. 3203*00 

19 

.98D4+Q0 

.3203+30 

. 3 2Q3 + OC 

.3203*30 

.3203*30 

.3203*00 

.3203*00 

18 

.9555+QO 

.3270+00 

. 3273+ 00 

. 3273 + 30 

. 3270+00 

.3270+00 

• 3270 + 00 

17 

.9237+qQ 

.3253*30 

.3253+00 

.3253+00 

. 3253+00 

.3253+00 

.3253+00 

16 

.8837+00 

.3226+00 

•3226+ GO 

.3225*30 

. 3226+03 

.3226+00 

. 3226*00 

15 

.8345+00 

.3180+00 

.3183+00 

.3180+00 

. 3180+30 

.3180+00 

.3183*00 

14 

.775b+0Q 

.31d9+0u 

.3109+00 

.3109+30 

. 3109+00 

.31 09+00 

. 3109+00 

13 

.7072*00 

.3032*00 

. 3C02+UO 

.3302*30 

. 3002+00 

.3302+00 

• 3002 + 00 

12 

.6307+00 

.2846+30 

,2646+00 

.2346*30 

. 2846+00 

.2846+00 

• 2846 + 00 

11 

.5487+00 

.2625+00 

.2625+00 

.2525+00 

.2625 +Qo 

.2&25+D0 

.2625*00 

IQ 

.4645+00 

.2324+00 

.2324+00 

.2324+30 

.2324+00 

.2324+00 

• 2324*00 

9 

• 3 8 1 6 + QO 

,1951+00 

.1931+00 

.1931+30 

. 1931+00 

.1931+00 

. 1931*00 

8 

.3040+00 

.1442+00 

, 1 442 + GC 

, 1442*00 

. 1442+00 

.1442*00 

• 1442*00 

7 

.2339+00 

-8 7o5-ol 

.8705-01 

.6705-31 

. 8705-01 

.3705-01 

.8705-01 

6 

,1730+00 

.2352-jl 

.2352-01 

.2352-01 

.2352-01 

.2352-01 

.?3^2-01 

5 

•1219+00 

.4237-01 

,4207-01 

.4207-11 

• 42n7-4U-_ 

.4207-01 

.423 7.:rJn 

4 

.8QDQ-01 

.1627+00 

.1627+00 

.1527*30 

. 1527+00 

.1527+00 

• 1627 + 00 

3 

.4655-01 

.3136+00 

.3136+00 

•3L36+3Q 

.3136. + a0.. 

. 3lX7* aD— 

.31 37*00 

2 

.2030-0 1 

.4092+02 

, 4092+ q2 

. 4392+32 

,4092+02 

.4092+02 

.4092+02 

1 

.3000 

"•927 2*02. 

-.9272+02 

. • 9?_7..2+_12._ 

9^.7 2 .+ Q 2 . .. 

-.9272+02 

-.9272+02 


o 

VO 



I 


H 

O 


JX = 

3, VARIABLE 

AT X ( 3 ) = 

.96821+00 




-- 




LEVEL 1 


tke 

♦ ♦♦ ♦ ♦ 



12 = 

1 

-> 

d_ 

3 

4 

5 

6 


Z - 

.3030 

. lOOO + QQ 

.2000+30 

• 3jgO+00 

.4000+00 

.BOOO+pQ 

20 

•9995+DQ 

.1531-j5 

..isai-Qi 

.1501-03 

. 1501-03 

.1501-03 

• 15nl-n3 

19 

.9&0M +D0 

.1631-33 

. 1501-03 

.1501-03 

.15 01-03 

.1501-03 

. 1501-03 

18 

• 9 b 5 5 + G 0 

.1531-03 

.1501-Q3 

.1501-03 

. 15 01-0 3 . 

.1501-03 

. ISDl -03 

17 

.9237+uq 

.1531-33 

. 15D1-03 

. 1531-03 

.15 01-03 

.15ol-03 

. 1501-03 

16 

.3837+00 

.1501-03 

.1501-03 

.i5oi-03 

.1501 -0 3 

. _ j I SJllr.D.I.- 

. 1501-03 

lb 

.8345+GO 

.16o1-q3 

.1501 -03 

. 1501-03 

.15 01-03 

.1501-03 

. i501-03 

1*1 

, 7756+00 

.1531-03 

. .1 5 01 r 0.3 

.1501-03 

. .I_5.air.a3_ 

.15Q1-03 

> 1SQ1-Q3 

13 

.7372 +00 

.1501-03 

. 1501- q3 

. 1501-03 

. 1501-03 

.1501-03 

. 1501-03 

12 

.&3C7+Qq 

.1501-03 

.1501-03 

.1531-33 

.1501-03 


. 1501-03 

11 

.6487+00 

.1501-03 

. 1501-03 

•1501-q3 

. 1501-03 

.1501-03 

. 1501-03 

ID 

.4645+00 

.1501-03 

. 1501-03 

. 1501-03 

. 1501-03 

.15nl-n3 

• 1501-03 

9 

.38 la + UQ 

.1501-03 

.1501-03 

.1501-03 

.15 01-03 

.1501-03 

. 1501-03 

8 

.3340+00 

.1500-03 

J-.^aQrOi 

.1500-03 

. 15 00-0 3 

.1500-03 

. 1500-n3 

7 

.,2 3 39 + 00 

• 1500-03 

. 1500-03 

. 1500-03 

. 1500-03 

.1500-03 

. 1500-03 

6 

. 1 7 30 + DD 

.1500-03 

. 1500-03 

.15J30-Q3 

. 1500-33 

.1500-03 

. 1500-03 

5 

.1219+GO 

, 1500-03 

. 1 5Q0-Q3 

.1500-33 

. 1500-03 

. 1 5qd-03 

. 1500-03 

M 

.Bono-oi 

.1500-03 

. 1 500-03 

.15QQ-33 

. l 5 .OD.ra_ 3 .. 

.1500-03 

. 1500-03 

3 

.4655-01 

.1484-03 

. 1 4S4-03 

.1484-03 

. 1484-03 

.1484-03 

.1484-03 

2 

,2o3o-Cl 

.2366-02 

.2383-02 

.2388-02 

.2388-n2 

.2388-02 

• 2388-02 

1 

.0300 

.2388-02 

. 23B8 -C2 

.2388-32 

. 2388-02 

.2388-02 

.2388-02 


(J 1 + ,^ 0 


JX- 

3, VARIABLL 

X 
1 1 

.96^^ Ol + OO 








L F V 8 L 1 


DIS 




IZ = 

1 

- 

3 

4 

5 

6 


1 - 

. 3 000 

• 1 L C.. 3 + u 0 

.2303+30 

• 3j3oa^30 

.4300+00 

.5000+00 

213 ■ 

.9995+00 

.671 7-0:4 

. 6 7 ] 7 - ,;,4 

.6717-04 

. 6717-04 

.5717-04 

.6717-04 

19 

.9804 +CU 

.6717-34 

. f 7! 7-04 

.6717-04 

.6717-04 

. *.6 7i7:-q4 

.6717-04 

16 

.9555+00 

.671 8-Q4 

71b-r4 

.6718-04 

.6718-34 

.5718-04 

.6718-04 

17 

• 9 2 3 7 + 0 D 

.6718-04 

,f-7ih-04 

.6718-04 

.6718-04 

.5718-04 

.6718-04 

16 

.8837+DD 

.6718-34 

» 6 / 1 b - 04 

.6718-34 

.67 l8 -04 

.5718-04 

.6718-04 

15 

.8345+00 

.6717-04 

. 6 71 7- ■:4 

.6717-34 

. 6717-04 

.5717-04 

.6717-04 

m 

. 7 7 5 b + c C 

.6717-04 

. e 717-04 

,u71 7-04 

.6717-04 

.5717-04 

.6717-04 

13 

.7070+00 

.671 6-J4 

. b71b- '4 

,67 i6-34 

.6716-04 

.5716-04 

.6716-04 

12 

.63C7+0D 

.6716-:,h 

• L llb~ ; 4 

.6715-34 

.6716-04 

.5716-Q4 

.6716-04 

11 

.5487+00 

. 6714-34 

.6714- ,4 

.D714-34 

.6714-04 

.5714-04 

.6714-04 

ID 

.4645+00 

.5713-04 

. 6 7 i 0 - . i 4 

.6713-34 

.6713 -04 

^ .6 713-_0_4_.__ 

.6713-04 

9 

.3813+30 


.6717-04 

. 6712-34 

.6712-04 

.5712-04 

.6712-04 

8 

.3040+00 

.6713-04 

.6710- JM 

.(^3 10-34 

.6710-04: 

..5J_l,D--a4_ 

_ .6710-04 

7 

.2339+00 

.67^9-04 

. 6704- C; 4 

.6739-34 

.6709-04 

.5709-04 

,6709-04 

b 

. 1 7 30 + Qu 

. b 7 0 P — ;J 4 

. r 7 Ob -04 

.67Q8-3 4__._ _ 

. 67 n8 -0 4 

.5708-04 

t6706-Q4 

5 

.1219+00 

» 6 7 3 8 - L.. 4 

, ^ 7 • P - 1,0 

• C70p-34 

.67QB-Q4 

.5708-04 

.6708-Q4 

4 

• 3 QDG~L) 

.6711-!’“ 

.6 7 10-;“ 

.6711-34 

.6711-24 

.b,7Ll-Jl4_ 

. f>7l 1 -04 

3 

.4S55“0i 

.6631-04 

. 9 L I - 3 4 

.0631-34 

.6601 -34 

.5501-04 

.6601-04 

2 

.2330 “■0 1 

, 4 4 4 [i - - . 

. 4 4 — O' ? 

. 4 ’4 4 3 - 3 2 

.4443-32 

. 4 44 n-n2 

.4440-02 

1 

0333 


. 6 f 7 9 - i 0 

.o879-q2 

,8879-32 

.8379-02 

.8879-02 



H* 

N3 


JXr 

3» VaRiAS 

LO at X( 5)- 

.96821+00 








LEVEL 1 


LrM 




iz- 

1 

2 

3 

4 

5 

6 


Z - 

.0000 

• 1 Q 0 0 + 0 Q 

._23Q0+_00. . 

^ODUtOQ 

.43QC+0D 

.5DQ0+QQ 

2D 

.9995+00 

.4500-02 

.*4 500^32 

,4500-02 

. 4500-02 

.4500-02 

.4500-02 

19 

•9BD4+DD 

. 453 Q -02 

.4500-02 

.4500-02 

.4500-02 

• 4 5 on-fl2 

. fttiDa-D2 

18 

. 9555 +DD 

.450D-Q2 

.4^00-02 

.4500-02 

. 4500-02 

.4500-02 

.4500-02 

17 

.9237+uC 

.4500-02 

.45PQ-G2 

.4500-02 

.4500-02 

.4500-02 

•45nn-o2 

16 

.88 37 + CO 

.4500-02 

.4500-02 

.4500-32 

•450q-P2 


. 45QD-D2 

15 

.83*45 + 00 

.4500-02 

.4500-02 

.4530-02 

.4500-0? 

.4500-02 

.4500-02 

1*4 

•7756+ UO 

.4502-02 

.4500-02 

.4500-02 

.4500-02 

•45oq-q2 

• 4500-02 

13 

.7072+00 

.4500-02 

.4500-02 

.45 Q0-n2 

. 45 C D — 0 2 

.4500-02 

• 4500-02 

12 

.6307+00 

.4500-02 

.4500-02 

.4500-32 

.4500-02 

.4500-02 

.4500-02 

11 

•5487+00 

.4500-02 

. 4 500-02 

.4500-02 

.4500-02 

.4500-02 

.4500-02 

10 

.46*45 + 00 

.4500-02 

.4500-02 

• 45 00”02 

.4500-02 

.4500-02 

. 4500-02 

9 

.3818+00 

.4500-02 

.4500-02 

.4500-02 

.4500-02 

.45on-D2 

.4500-02 

8 

.3Q4Q+D0 

.4500-02 

.4500-02 

.4500-02 

• 45 00-0 2 

,4500-02 

.4500-02 

7 

.2339+00 

.4530-02 

.4500-02 

.4500-02 

.4500-02 

.4500-02 

.4500-02 

6 

. l7 30 + CO 

.4500-02 

.4500-02 

.4500-02 

•4500-02 

.4500-02 

.4500-02 

5 

.1219+00 

,4530-02 

.4 5 00-02' 

.4500-02 

. 4500-02 

.4500-02 

.4500-02 

*4 

• 8 0 00“Gl 

.4500-02 

•4500-02 

.4500-02 

.4500-02 

.4500-02 

.4500-02 

3 

.4655“ul 

.4499-02 

. 4499-02 

.4499-02 

.4499-02 

.4499-02 

.4499-02 

2 

.2030-01 

,4318-02 

.4318-02 

.4318-02 

.4318-02 

.4318-02 

.4318-02 

1 

• 0 0 G 0 

.0000 

.0000 

.0000 

. DODD 

.0000 

.0000 



113 




JX- 

3t Variable 

AT X ( 3 )- 

. 968 21 + 00 


- — 

— 




' ' ■ 

level’" 1 ' 

4c 4c 4: 4c 4; 

T AU 





1 


3 

4 , 

5 

6 


z - 

. J D 0 0 

. ! 200 + OG 

.2300+00 

. 330C+00 

.4000+00 

. 5000+00 

2D 

.9995+DD 

.4047-01 

,4047-01 

.4047-01 

. 4047-01 

,4047-Dl 

.4047-01 

19 

.98G4+LD 

,452:0-01 

.4520-01 

.4520-01 

.4520-01 

•4520-01 

.4520-01 

16 

•9555+UC 

.4614-01 

.4214-01 

.4514-01 

.4614-31 

.4514-01 

.4614-01 

17 

.9237 + IjO 

,4590-01 

.4590-01 

.4593-01 

.4590-01 

.4590-01 

.4590-01 

16 

.B837+0C 

.4551-01 

.4551-01 

.4551-01 

.4551-31 

.4551-01 

, 4551-01 

15 

.3345+00 

.4486-01 

.44SS-Q1 

.4486-11 

.4,18 6 

-44ft6-ni 

■44R6-0I 

14 

.7756+00 

,4386-01 

.438b-cl 

.4386-31 

.4386-31 

.4386-01 

.4386-01 

13 

•7D72+uO 

.4236-01 

,4236-01 

•4236-31 

.42 36.-11 . 

■ ■■,.•■4 2 3 6^.^1 

.4236-^1 

12 

•6307+00 

.4015-01 

.4015-01 

.4315-01 

.4q 15-31 

.4015-01 

.4015-01 

11 

.5487+00 

,.37Q4-ai 

. 3704-01 

,..3I 

,x37QiL-Q-I_ 

■■.3 7o4-gi 

. 3JLD-4-01 

ic 

, 4 6 4 5 + Q C 

.3279-01 

.3279-01 

.5?79_01 

. 3279-01 

•3279-01 

•3279-01 

9 

.3818+00 

.2724-01 

.2724-01 

.2724-01 

.2724-3 1 , 

»?7?4-ni 

■ 2724-01 

6 

.5q4u+0D 

.2C34-01 

.2034-01 

.2334-31 

. 2034 -31 

.2334-01 

.2034-01 

7 

.2339+GD 

.1228-01 

.1228-01 

* 1228-01 

. 12 2E-11 



6 

•173Q+00 

.3319-02 

.3319-02 

, 3319-02 

. 33 19-Q2 

.3319-02 

.3319-02 

5 

.1219+00 

.5936-02 

."936-02 . 

..S?16-J32... 

. 59 3.6_-Q.^.. 

..5936-JI2 

.5936-Q2 

4 

.8DDQ-j;;;i 

.2295-01 

.2295-01 

.2295-01 

.22 95_31 

.2295-01 

.2295-01 

3 

.4655-31 

.4428-01 

.4423-01 

.4428-31 

. 4428-01. 

._4 4.2B.-d]! 

_^4jL28-Dl 

2 

.2033-01 

,6129+01 

.6129+01 

.6129+01 

.6129+01 

.5129+01 

.6129+01 

1 

.0330 

.2411+02 

.2411+02 

.2411+32 

.24ll-»02 

.2 4.11+02 

.2411+02 



M 

I-* 


JX- 3* variable at X( Z)- . 9682,1+00 





LEVEL 

] ^ 4c t 4c 

CP/2 




12- 

1 

O 

c 

3 

4 

5 

6 


Z : 

.3000 

. 1033+ DC 

.2330+30 

* ^DOD+OO 

.43QC+0Q 

•50QQ*0Q 

20 

.9995+00 

.2 193-01^ 

.2217-06 

.21 98-36 

.2105-0 6 

.2272-06 

.2272-06 

19 

.9804+00 

- .7916-36 

-.7991-06 

-.8384-36 

83 p6 -o6 

- . 7g2 3-06 

-.7935-06 

18 

•9555+ QG 

- .8939-05 

-. 8894-05 

-.8348-35 

-.8866-05 

-.8883-05 

-.8883^05 

17 

.9237+00 

-.3466-35 

-.8464-05 

-.8451-05 

-.3533-05 

-.8494-05 

-. 8494-05 

16 

.8837+00 

-.7970-35 

-.7989-05 

-.Bo ‘*5-3.5 

-.£317-05 

-.3300-05 

-. 8034-05 

15 

.8345+00 

-.6979-05 

-. 6991 -05 

-.7317-35 

-.7121-05 

-.6983-05 

-.■^002-05 

14 

,7756+00 

-.6398-05 

-.6419-05 

-.6436-05 

-.6448-n5 

-.6402-05 

-.6426-05 

13 

.7072+00 

-.6270-05 

-. 6275-.05 

-.6?9e-05 

-.6296-05 

-.6296-05 

-.6298-05 

12 

.6307+00 

-.6586 —'0 5 

-.6586—05 

-.6531-35 

-.6545-05 

-_6545-a5 

-.■654,2 -^5 

11 

.5487+00 

-.6836-05 

6914-05 

-.7300-05 

-.6920-05 

-.6979-05 

-.6963-05 

10 

.4645+00 

- ,39i 3-Utj 

-.8937-05 

-.9313-05 

-.8937-05 

-«9nn8-n5 

-.9008-05 

9 

.3318+00 

-.1160-34 

-. 1159-34 

-.1156-34 

-.1160-04 

-.1157-04 

1156-04 

8 

.3040+00 

“ • 1 57 o~04 

-. 1566-04 

-.1564-34 

-. 1562-04 

-.1 5 64-n** 

156c;-04 

7 

.,2339 + DD 

-.2077-04 

-. 2074-04 

-.2366-34 

-.2d73-34 

-.2365-04 

-.2066-04 

6 

•173^+00 

-.2468-34 

-,2466-04 

-.2460-34 

-. 2464 -34 

-.2466-04 

-.2464-04 

5 

.1219+ 00 

-.3115-04 

-.3112-34 

-. 51 03-04 

-.31 09^0 4 

-.3107-04 

-. 3107-04 

4 

.630u-U1 

,8065-04 

.8065-04 

,8362-34 

. 8D67-D4 

,8366-04 

.8065-04 

3 

.4655-01 

.1052-03 

,1052-03 

.1351-03 

. 1D51-D3 

.l051~03 

. 1051-03 

2 

.2 j3Li-Dl 

,391 3-32 

,3913-02 

. 3913-02 

.3913-02 

.3913-02 

.3913-02 

1 

.0000 

,5576-02 

.5575-02 

.5575-02 

.5575-02 

.5575-02 

,5575-02 





JXr 

3, VARIABLE 

AT X ( 3 ; 

f- .96^21 + 00 









LEVEL 1 


MA CH _ 




IZ = 

1 

2 

3 

4 

*5 

6 


2 = 

.3300 

. 1'003+GC 

.2030+00 

. 33G0+QG 

.40Q0+D0 

.5000+00 

2d 

,9995+CJ 

.3001+31 

. 3C01 + G1 

.5001+01 

• 3DQi+ni 

iiLlQl+_0J 

.3001+01 

19 

.98C4+U0 

.J'jJl + 31 

.7031+01 

.5001+31 

.3001+01 

.3001+31 

.3001+01 

18 

• 95 55 + CG 

.3Qj 1 + Qi 
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